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PREFACE 

Reliability models generally deal with two types of problems depending 
upon system stages viz., design and operative stage. Problems in design stage 
pertain to redundancy allocation, reliability maximization and/or cost 
minimization etc., whereas problems in operative stage deal with analysis of 
stochastic behaviour of the systems and to determine optimal maintenance 
policies. 

The review of existing literature on reliability reveals that authors who 
have contributed to the economic aspects of reliability had definite cost 
structures in their minds to superimpose on the system model. Generally 
speaking, reliability and cost are related parameters of a system and in fact 
other parameters of the system can be expressed in terms of cost. 

It is needless to mention that problems at both these stages are 
important in their own right and require due attention. However, for the 
analysis of maintained systems, the problems at the latter stage (operative) 
play a dominant role because a system analyst loses his control over the 
design parameters once the system enters operative stage. So, for maintained 
systems there is a genuine demand of research investigation in the direction 
of repair/maintenance which can help a maintenance engineer/analyst in making 
an optimal choice of repair rule (for a fixed failure rate) and/or in having 
better understanding of system’s economics by studying parameters’ effectiveness 
from cost considerations. 

The thesis consists of four chapters ; 
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Chapter I : Introduction 

Chapter II : Optimal Repair Rate in Some Maintained Systems 

Chapter III : Optimal Failure Analysis in Series Systems 

Chapter IV : Evaluation of Series System Through Dynamic Programming 
Technique. 

Chapter I presents a critical review of the existing literature and a 
brief account of contribution made by the author. 

Chapter II discusses the profit evaluation studies for some systems with 
repair maintenance. Three models describe single unit maintained systems 
under different operational conditions and one model deals with two-unit 
system with parallel redundancy. The concept of minor repair has also been 
applied in one of the models. The cost structure of Howard has been 
superimposed on the Markov process / semi-Markov Process generated by a 
system model. The expected profit thus obtained has been maximized and 
optimal repair rates are determined. This cost structure allows different cost 
rates (variables) in each state and different transition costs incurred at the 
time when the system makes a transition from one state to another. Following 
4 models in this chapter are discussed in the following sections : 

2.2 Profit evaluation in a single unit system 

2.3 Profit evaluation in a single unit system with a less productive 

state. 

2.4 Profit considerations in a single unit system with minor repair. 

2.5 Profit evaluation in a 2-unit parallel redundant system. 



Chapter III deals with failure analysis in series systems. Problems of 
failure classification and resource allocation on failure detection mechanisms 
have been considered. For the problem of failure detection in a series system, 
a cost structure involving probabilities of correct detection (depending upon 
cost spent on detection) has been mixed with Howard’s reward structure and 
is superimposed on a 2-unit series system. Optimal allocation of cost on 
detection mechanisms has been determined when the total cost for detection 
is fixed. Further a problem of failure classification in a series system has been 
discussed after taking into account cost of misclassification whenever the system 
fails. A decision rule based on the life-time distribution of the system which 
minimizes the expected misclassification cost has been suggested. A procedure 
has also been suggested to examine the discriminatoiy power of the rule. 
There are following two formulations discussed in this chapter : 

3.2 Optimal failure classification in an n-unit series system. 

3.3 Optimal allocation of cost to detectors in a 2-unit series system. 

Chapter IV discusses application of dynamic programming technique to 
solve problems in series systems. Problem of cost allocation has been formulated 
and solved through dynamic programming technique. This is discussed in the 
following section : 

4.2 Optimum cost allocation in a series system under a constraint 
on the cost. 

It may not be out of place to mention that the results contained in 
different chapters are based on seven research papers published in different 
international journals of repute by the author. It may be mentioned that the 
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procedure followed for identifying the equations in the thesis is to put numerals 
within parentheses, where the first figure refers to the chapter, the second 
represents the section number and the third represents the serial number of 
the equation, e.g., (3.2.5) denotes the 5th equation of section 2 of chapter 
III. The same pattern is followed to number the tables and figures also. 

Abbreviations : 

s-independent : statistically independent 
t : monotonically increasing 
I : monotonically decreasing. 
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CHAPTER I 

INTRODUCTION 



1.1 INTRODUCTION 


The reliable performance of system is of immense interest in modern 
business, industrial and defence systems. In fact, increasing complexity of 
present day system has resulted in engineering problems involving high 
performance, reliability and maintainability. It is the importance of reliability 
in system analysis studies that has attracted attention of a large number of 
research workers from different disciplines to cany out their analysis from 
the reliability view point. Investigation of reliability models with cost 
considerations is of utmost importance to system analysts. 

We develop below basic definitions and terminology which are needed 
to understand a reliability model. 

1.2 Definitions and Concepts 

The terms/parameters which are used in reliability may be divided into 
two categories : 

(A) System configuration parameters 

(B) System effectiveness parameters. 

(A) System Configuration Parameters : 

1. Series system is one in which a failure in any one of its units 
(components) leads to system failure. 

2. Redundant systems are those where more than one means 
(units) are available for performing a function, so that when one 
means (unit) fails, others are available to perform the task. 
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Parallel or active redundant sjstems are those where all 
redundant units and basic units operate simultaneously rather 
than being switched on when needed. 

4. Standlty redundant systems are those in which only one unit 
operates at a time, while the other redundant units (components) 
are switched on sequentially only when needed i.e., upon the 
failure of the operating units. Here, the system is said to have 
failed when all the units have failed. 

Tlie standby redundancy may be cold, warm or hot. 

(a) Cold standby redundant systems are those in which redundant 
units do not fail when they are in standby. 

(b) Warm standby redundant systems are those where redundant 
units can fail when they are in standby. Normally, in such systems failure 
probability of standby unit is less than that of an operating unit. 

(c) Hot standby redundant systems are those in which redundant 
units can fail when they are in standby with failure probability same as that 
of an operating unit. 

5. k-out-of-n : F system it consists of n units and the system 
fails if and only if k (< n) out of its n units are failed. Such 
systems are prevalent in practical systems, e.g., a communication 
system with 4 transmitters in which the average message load is 
such that at least 3 transmitters must be operational for successful 
message delivery, thus the system is 2-out-of-4 : F system. 
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6. k-out-of-n : G system it is a redundant system composed of 
n units and the system functions if and only if at least k ( < n) 
out of n units are operable. It is clear that series and parallel 
systems are n-out-of-n : G and 1-out-of-n : G systems respectively. 

7. Complex systems are those systems which do not fall into the 
category of either standby or parallel or series-parallel systems. 

8. Imperfect switchover The phenomenon which occurs when a 

standby unit is switched to operate at the time of failure of an 
operative unit in a redundant system is referred to as ‘switchover’. 
Some times the switchover is not perfect i.e., not 100% error 
free and/or instantaneous. This is called ‘imperfect switchover’. 

9. Maintainability may be defined as a characteristic of design 
and installation which is expressed as the probability that a unit 
will be retained in or restored to a specified condition within a 
given period of time, when the maintenance is performed in 
accordance with certain rules. Maintenance includes all actions 
necessary for retaining a unit in or restoring it to a specified 
condition. 

It is fiirther divided into following two categories 

(a) Repair maintenance whenever a unit fails, an activity is performed 
to restore the unit to a specified condition (operation). 

(b) Preventive Maintenance (P.M.) this is performed in an attempt 
to retain a unit in a specified condition by providing systematic, inspection. 
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detection and prevention of failure. 

P.M. is required for costly systems like computers, satellites, radar, 
transmitters etc., where cost of failure or replacement is quite heavy and 
systems are required to operate without failure for a long time span. There 
are several types of P.M. policies. 

10. Service replacements are those replacement which are made 
after the occurrence of failure. 

11. Planned replacements are those replacement which are made 
after the occurrence of failure. Such replacement may be made 
at scheduled times. 

12. Delayed repair the repair of a unit is some times delayed due 
to any physical or environmental reason. This process is called 
delayed repair. 

13. Fault detector occurrence of a fault (or failure) in a system 
is not always detected automatically. As such some mechanism 
may be used to detect the fault, this is referred to as fault 
detector. 

(B) System effectiveness parameters : 

1. Reliability may be defined as the probability that a system 
will perform its function adequately for a definite period of time 
under specified conditions. 

2. Pointwise availability probability that the system will be able 
to operate within certain tolerances at a given instant of time. 
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3. Interval reliability probability that a system continues to 
operate during a given interval of time. 

4. Interval availability expected proportion of a given interval 
of time that the system will be able to operate within certain 
tolerances. 


5. Steady-state availability expected proportion of time that a 
system is operable when it is allowed to operate for a very long 
period. 

6. First passage time to system failure this is defined as the 
probability that starting in a given state (up) system enters a 
failed state for the first time in time t. 

That is 

^of( 0 = P[system goes to state ‘F’ in time < t [system is in ‘O’ at t = 0]. 

7. Mean time to system failure if T denotes time -to -first -system 
failure starting from the beginning, then MTSF is given by E(T) 
where E stands for mathematical expectation. 

This measure of system effectiveness appears to be more appropriate 
for the so called one-shot systems, e.g., missiles, rockets etc., which are used 
once for all. 

Infact, first passage time distribution to system failure, ^^(t) and MTSF 
are related to the following manner. 


MTSF = 


5=0 
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where (pl{s) denotes L.S. transform of ^o(t). 

8. Expected profit This is the total expected earnings per unit 
time of the system if the system is allowed to operate for 
sufficiently long period. 

13 System Models in the Field of Reliability 

Keeping in view the varied nature of practical problems, a large number 
of reliability models have been developed and analyzed in literature. 

Most of the papers appeared are motivated to the following two 
directions:- 

(1) Analysis of special systems. 

(2) Development of reliability, availability and cost optimization models. 

In the analysis of special systems one is interested in a specific system. 
The general approach for the analysis of such systems consists in defining and 
formulating the system model i.e., to state its assumptions, to enumerate 
different states and identify them as up and down ones, to develop the 
underlying equations governing the system behaviour and to apply an appropriate 
technique to obtain the solution. 

For the second case viz., in the development of reliability, availability, 
cost optimization models, first of all one has to formulate the problem clearly 
i.e., to construct the objective function and constraints. Infact, in this area, 
formulation aspect is important in its own right, this involves clear understanding 
of effectiveness criterion, constraints visualisation, and their impact on objective 
functions. Once the problem is formulated successfully, next step is to search 
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and/or develop a suitable solution method to get optimal values for decision 
variables. 

We now review the existing literature on reliability keeping in view 
the above classification. 

(1) Analysis of special system In general following types of system 
models have been discussed under this category. 

(i) Cold standby redundant systems Muth (1966) cosidered a 
cold standby system with constant hazard rate for the on line 
unit where the repair time is a random variable having a gamma 
distribution and obtained expressions for the Laplace transform 
of mean time-to-system-failure (MTSF) of the system. 
Subsequently, Linton and Baraswell (1973) have generalized the 
results of Muth (1966) for case in which life time of the unit 
and its repair time are generally distributed random variables. 
Nakagawa and Osaki (1974a) obtained for the same system, 
expressions for the Laplace transform of the reliability, steady-state 
availability and its MTSF. Nakagawa and Osaki (1974b) derived 
expressions for the expected number of visits to a state. Srinivasan 
and Gopalan (1973a) obtained expressions for the Laplace 
transform of reliability and pointwise availability. Kumar (1977a) 
obtained expression for steady-state profit for a 2-unit cold 
standby redundant system. 

(ii) Warm standby redundant systems :- Branson and Shah (1971) 
considered a 2-umt system with exponential failure times and 
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general repair times and obtained the MTSF and steady-state 
availability. SriraVasan and Gopalan (1973b) obtained the Laplace 
transform of pointwise availability of the system comprising of 
dissimilar units by employing supplementary variable technique. 
Kumar (1977d) considered a 2-unit system with exponential 
failure times and general repair times and obtained the expression 
for expected profit of the system. Birolini (1975) analyzed his 
model when the failure time distribution of the on line unit 
follows general distribution and failure time of the standby unit 
follows exponential distribution, repair times are also of 
exponential nature by employing regenerative stochastic processes. 

(iii) Parallel redundant systems :- Gaver (1963) studied a 2-unit 
parallel redundant system with constant state dependent hazard 
rate and arbitrary repair and obtained the Laplace transform of 
reliability of the system its MTSF and steady-state unavailability, 
using supplement variable technique. Kulshrestha (1968, 1970), 
Mine and Kawai (1974), Linton and Saw (1974) have studied 
some parallel redundant systems. Kodama and Deguchi (1974) 
obtained the MTSF of a parallel redundant system with Erlangian 
failure and general repair. Linton (1976) extended these results 
for the case when the failure and repair distributions of one of 
the units are Erlangian and those of the other are arbitrary by 
employing supplementary variable technique. 

(iv) 2-Unit standby redundant system with imperfect switch over :- 
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The importance of switch in standby redundant systems has 
attracted the attention of many research workers to incofporate 
switch behaviour in investigating the stochastic behaviour of 
standby redundant systems. Gopalan (1975), Osaki (1972a) 
Prakash (1973), Nakagawa and Osaki (1975), Nielsen and Runge 
(1974), Kumar (1977b) have studied such systems under different 
operational circumstances. Chow (1972, 1973) constructed 

mathematical models for the reliability of modularly redundant 
systems with unequal failure rates for the operating and standby 
units and different types of svdtch failures. Arnold (1973) 
introduced the concept of ‘coverage’ for repairable systems. 

(v) Redundant systems with delayed repair When a unit or switch 
fails; repair facility may not be available because of its 
preoccupation with higher priority jobs and as such failed unit 
might have to wait for sometime. Several other reasons may also 
be attributed to delay in repair. The reliability models with delay 
in repair have been analyzed by a large number of workers. 
[Gopalan and Saxena (1977), Kapur and Kapoor (1975a), Khalil 
and Bogus (1975), Khalil (1977b), Kumar and Lai (1979) etc.]. 

(vi) Redundant systems with two types of repair Some electrical 
and mechanical devices can fail in two or more mutually exclusive 
failure modes. A vending machine or public telephone fails due 
to failure of units under benign environments or to the failure 
of human misuse. The failure (and thus equivalently repair) may 
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be categorized in 2 classes.Many papers with 2 types of repair 
have appeared in the literature, [Garg and Kumar (1977), Khalil 
(1977a, b), Kumar and Agarwal (1978), Osaki and Okumoto 
(1977), Proctor and Singh (1976) etc.]. 

(vii) Complex ^sterns All redundant systems may not be either 
standby or parallel or series-parallel. Consider a complex system 
consisting of two classes of components denoted by Lj and 
Class Lj has n components in series and class Lj identical 

components in parallel redundancy. Such systems have been 
investigated by Garg (1963a,b), Gupta (1973a,b), Sharraa (1974), 
Agarwal (1975). Such systems have also been considered by 
Varma (1973) with repair priority and by Kapur and Kapoor 
(1975b) with demand pattern. 

Different techniques were used by research workers in the evaluation 
study of reliability models to obtain parameters of interest. For example 
Gnedenko et al. (1965), Osaki (1970b) and Buzacott (1971) have used the 
renewal theoretic approach for the analysis of their models. The technique 
of Markov-renewal process has been successfully employed for the analysis of 
redundant systems by Srinivasan (1968), Osaki (1970a, 1972a,b), Osaki and 
Okumoto (1977). Further, Nakagawa and Osaki (1974a) have suggested a 
unique modification of MRP which is useful for the analysis of 2-unit redundant 
systems under general failure and repair time distributions. Later, this technique 
was employed in a large number of papers. [Nakagawa and Osaki (1974b, 
1975), Kapur and Kapoor (1975a,b), Kapil and Sinha (1981)]. Gnedenko (1967), 
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Srinivasan et al. (1973) have used regeneration point technique for the analysis 
of redundant systems. The supplementary variable technique has been used 
in standby redundant system by Gaver (1963) and subsequently by others, 
Garg (1963a,b), Garg and Kumar (1977) to treat general distributions for 
repair times. Various other methods have also been applied by workers to 
analyze their models, e.g., discrete transforms Kulshrestha (1972), conditional 
transforms, Linton and Braswell (1973). 

Branson and Shah (1971) applied the results from the theory of 
semi-Markov processes to reliability analysis for a 2-unit warm standby systems. 
Later on, this was applied to other systems [Kumar (1976, 1977d), Kumar 
and Lai (1979,1980)]. 

But a close look at the processes considered in these papers reveals 
the fact that the underlying process is not SMP in true sense. However, it 
looks so. The interesting fact is that the process satisfies certain conditions 
[Arndt (1977), Arndt and Franken (1977)] under which results from SMP 
theory are applicable. This has been pointed out in Kumar et al. (1981). 

(2) Development of Reliabilitity, Availability and Cost Optimization 
Models :- A large number of developments have appeared to optimize system 
effectiveness: reliability, availability, cost, profit etc. For problems of this 
nature, Tillman et al. (1980) may serve as a good reference. We review below 
some of the important optimization models. 

(i) Optimization models in series and parallel system In these 

models there are several units (stages) say, connected in series 
and the units are statistically independent. If reliability of the 
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unit is denoted by Rj, the system reliability R is given by 

R = n Rj 

j=i ^ 

Many workers [Ghare and Taylor (1969), Beraha and Misra (1974)] 
have made attempts to develop algorithms which maximize R subject to k 
resource constraints of the form 

i hij(Rj) < a. i = 1, 2, k 

j=i 

where hjj(Rj) denotes the requirement of the i*'’ resource and a-^ stands for 
availability of the i'*’ resource. The functions hij(Rj) may be linear or nonlinear 
depending upon situations under investigation. 

Another useful area in series optimization models is redundancy 
allocation in an optimal manner. Banerjee and Rajamani (1973), Shetty and 
Sengupta (1975), Luus (1975) have considered problems to choose the optimum 
number of components, say mj in parallel at each of the n stages in a series 
system so that reliability is maximized. Mathematically, the problem in general 
form is 

Max R ; R = n (1 - QfO, Qj =1-Rj 

j=i 

subject to the constraints given earlier. 

Misra and Ljubojevic (1973), Tillman et al. (1977b) considered a 
simultaneous reliability and redundancy allocation problem viz., 

Max R ; R = n [1 - (1 - 

j=i 
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subject to 

n n 

1 hij(Rp = 2 hjj(R.,m.) < 

j=i j=i 

and choose R^ and nij. 

They have used Lagrangian Multiplier technique to solve the problem. 
Shetty and Sengupta (1975) applied slacked sequential unconstraint minimization 
technique (SLUMT) which is superior to SUMT used by Tillman et al. (1970) 
for solving their problems. 

Another formulation in series models has been due to Aggarwal (1977). 
He discussed the problem of minimizing the total cost under a lower bound 
on system reliability. In other words, 

n 

Min C ; C = X C/mj) 

j=i 

subject to 

R = fi Rj(mp > Ro 

j=i 

Cj(mj) being the cost of the j**" stage when m^ parallel components are 
used at this stage. He used heuristic approach for the solution. 

Misra (1971a) gave a dynamic programming formulation of redundancy 
allocation at each stage under two nonlinear constraints. 

The problem of choosing MTBF, M l TR and mj by cost minimization 
subject to availability constraints has been formulated by Lambert et al. (1971). 
They have taken total cost C and availability A for a parallel series system 
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as 


C = i (R,,M,,mj) 

j-i 

n 

A = n Aj(Rj,Mj,mj) 

A --5^ 

^ Rj+Mj 

Rj and Mj are MTBF and M l I K at the stage respectively. 

The problem is to minimize C by choosing Rj, Mj, m^ within permissible 
design requirements so as to meet a specified level of availability. 

A typical design requirement for the stage is 

R. > 4mj , Mj < 5 

Dynamic programming technique was used to obtain the solution. 

(ii) Reliability optimization models with repair schemes Pal and 
Bhattacharjee (1978) considered a system having n stages and 
stage has mj units. The objective is to maximize system 
reliability subject to nonlinear constraints depending upon 
maintenance. The system reliability is 

R(t) = n Rj(t) = exp (-Zt) 

j=i 

where 

z = i z 

j=l 
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and 


l-/)p[l+m.(l-pp] 


- failure rate of stage j _ ^ 

repair rate of stage j 

The steady-state probability Pqj of stage j with its units being idle is 



The problem is 
Max R(t) 

m. 

J 


subject to 

mj > 1 (j = 1. 2, n) 

n 

2 hjj mj < Cj (resource constraints) 

j=i 

t(l-Poj) T*” (busy period constraints) 

(iii) Optimization in repair, replacement and inspection models 

Replacement is carried out when the failures cause serious 
damages and/or imply expensive repairs. Replacement can be 
made in accordance with a predefined policy. 

Replacement theory is of great importance in reliability theory and has 
been investigated by many research workers. Barlow and Proschan (1965) 
discussed several interesting replacement problems and obtained the optimum 
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repair policies minimizing the expected cost per unit time. Glasser (1967) 
considered the age replacement problem of Barlow and Proschan (1965) and 
gave graphs of the optimum repair policies. Fox (1966) introduced discount 
rate for a replacement problem. Scheaffer (1971) investigated a replacement 
model with an increasing cost for an operating unit. Nakagawa and Osaki 
(1974c) obtained the optimum policies for a replacement model with delay. 
Nakagawa and Osaki (1974d) have taken a repair limit replacement policy 
which is defined as the policy in which the repair is stopped if it is not 
completed within a fixed time (called the repair limit time). The problem is 
to determine the optimum repair time limit r, which minimizes the expected 
cost per unit time for an infinite time span. Kapil and Sinha (1981) have 
discussed repair limit suspension policy for a 2-unit redundant system with 
2-phase repair. A failed unit is first repaired by a phase 1 repair. 

If the repair is not completed in a fixed time then the unit goes for 
a phase 2 repair. They have obtained optimum repair limit suspension policies 
which maximize the availability using MRP results. 

Muth (1977) also studied a similar problem in which a replacement 
takes place at the first failure time after a predetermined age t. 

Mine and Kawai (1975) discussed an optimal inspection and replacement 
policy for a unit which assumes any one of several Markov states. The policy 
evaluation function is expected cost per unit time over an infinite time span. 
Hie problem is formulated as a SMDP with a modified policy improvement 
routine. 

Alam and Sarma (1974) considered maintenance policies for a machine 
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with degradation in performance with age and subject to failure. 

Hastings (1969) developed a repair limit replacement method. In this, 
when an item requires repair it is first inspected and the repair cost is 
estimated. If the estimated cost exceeds a certain amount known as the repair 
limit, then the item is not repaired but replaced. Two types of problems, one 
in which condition of equipment is related to age and other in which condition 
is related to number of major repairs an item has had, have been analyzed. 
Dynamic Programming methods have been used to obtain optimum repair 
limits. 

Diveroli (1974) has considered simultaneously the problems of preventive 
replacement and service replacement. Several procedures have been developed 
to determine optimal policies. He also proves the optimal preventive 
replacement policy is a function of the optimal repair limit. 

(iv) Profit optimization in redundant system models Expected 
steady-state profit in redundant systems has been obtained in 
several papers [Kumar (1976), Kumar and Lai (1979)]. This 
measure of system performance was used to determine optimal 
P.M. in Kumar (1976). However, profit optimization problem 
was considered in Kumar and Lai (1980) to determine optimal 
maintenance policy for a redundant system. This method is based 
on Howard’s (1971) policy iteration. Srinivasan et al. (1971) have 
considered expected profit in a 2-umt redundant system by mixing 
of two renewal processes. Further, Kumar and Kapoor (1981a) 
extended the solution method to find optimal maintenance policy 


for a redundant system when the cost structure is discounted 
one. 

(v) Optimization models with fault detection The concept of fault 
detection in series system is important and needs attention. Gross 
(1970) had studied a 2-unit series system with fault detector, 
misclassification cost. He had minimized the cost of 
misclassification and determined the optimal cost spent on 
detector. Later, Kumar (1975) had extended the result to the 
case of n-unit series system. He obtained a sufficient condition 
under which a detection mechanism would be economically 
feasible i.e., expected cost of misclassification is minimized. 
Recently, Takami et al. (1978) have developed a Markov model 
for a class of series systems which have fault detectors to find 
component failures. The optimal allocation of fault detectors is 
determined by formulating the problem as 0-1 integer 
programming problem. 

Recently, Kumar and Kapoor (1981b) have solved a problem concerning 
classification of failure in a 2-unit series system, the classification criterion 
depends upon the time to system failure. 

1.4 Some Cost Structures in Reliability Models 

The review of existing literature in previous sections reveals that authors 
who have contributed to the economic aspects of reliability had definite cost 
structures in their minds to superimpose on the system model. Generally 
speaking, reliability and cost are related parameters of a system. But a unique 
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expression of this relationship which is valid in reliability modelling is not 
feasible mathematically. Never the less, some properties originating from the 
conceptual frameworks which are desired of cost fimctions are 

(I) Components with very low reliability cost very low. 

(II) Similarly, components which are highly reliable cost very high. 

(III) Cost is an increasing function of reliability. 

(IV) Derivative of cost with respect to reliability is an increasing 
function of reliability. 

We now discuss some important cost structures which have appeared 
in reliability modelling and serve a useful purpose to system analysts 

(1) Beriphol (1961) has expressed the cost of a component as a function 
of its reliability as given below 

C„ = YZ^^exp[-K,„(l-R„(C„))] 

where 

C„ : Cost of the n*** component 

R„(Cn) : reliability of the n*’’ component 

Kj„,K 2 n : constants depending upon components. 

A modified version of the above cost structure has appeared in Aggarwal 
and Gupta (1975) 

c. = K. [tanCfR.)]"^' 
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here f(Rn) is a function of R^. In many situations f(Rn) takes the form 


f(R„) = 1 + R.'. . 0 < X. < 1 

f(Rn) = . 0 < ni„ < 2 

(2) In the problem of maximizing R defined on page 12 

R = n [l-a-Rfi] 

j=i 

subject to 

i ^j(Rj) < a, 

j=i 

and to choose mj and Rj simultaneously, Misra and Ljubojevic (1973) have 
taken the functions hij(Rj)’s as 

hjjCRj, mp = a^j exp[-bi/(l - R^)] 

Infact most of the models assume additivity of costs i.e., total cost C 
is additive in terms of the costs of the constitutent units 

C = i q(m.) 

j=l 

where Cj(mj) is the cost of the j* stage assumed to be a known function of 
mj- 

(3) Many developments are based on expected cost per unit time for 
an infinite time span. Nakagawa and Osaki (1974d), Kontoleon (1977), Muth 
(1977), Mine and Kawai (1974, 1975) have used this concept in their works. 
This is obtained as follows. 
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Find the expected cost from the instant a system starts operating upto 
the completion of the repair/replacement. (This constitutes one complete cycle). 
Divide this quantity the mean time of one cycle. 

For example, if F(t) is cdf of repair time of a system, Cq is the 
replacement cost, Cj.(t) is the expected repair cost incurred during (0,7] then 
the expected total cost over one cycle is 

c = [q, + C,(r)] F(t) + / C,(t)dF(t) 

0 

= C„F(t) + 1 F(t)dC,(t) 

0 

where r is repair limit time. 

TTie mean time of the cycle is 

A+TF(t)+} t dF(t) = A+/ F(t)dt 
0 0 

Thus the expected cost per unit time over an infinite time span is 

CoF(t)+/ F(t) dC,(t) 

Ur) = f= 

A+/ F(t)dt 
0 

For Cj.(t) the two usual forms are : 

C,(t) = at^ a > 0, 0 < b < 2 
i.e., repair cost is proportional to time and 
C,(t) = a( e*** -1), a >0, b >0 

i.e., exponential repair cost. 
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The problem is to obtain optimal r which mim’mizes L(t) 

(4) Shershin (1970) has discussed a problem of simultaneous 
apportionment of reliability and maintainability. He visualized the total cost 
function for a system with N subsystems as the sum of the failure rate 
and maintainability costs. Assuming exponential failure rates, the failure rate 
Aj of the system and the failure rate for the subsystem are related by 

.N 

A. = 2 A, 

k=l 

Further, let d^(A^) be the cost for achieving failure rate for the 
subsystem. Then, the cost of achieving a system failure rate A^ is 

C^,.(AJ = 2 d,(A,) 

k=l 

He has concentrated on maintenance time and has considered 3 major 
activities and their associated costs to be performed under maintenance time. 

(i) Maintainability during design stage; 

(ii) Performing maintenance to prevent failures termed as P.M.; and 

(iii) Repairing a failed subsystem/subsystems termed as C.M. associated 
with each component. 

(5) Nakagawa and Osaki (1976a) considered a problem of optimization 
for a one-unit reliability model which is provided with n spares. They assumed 
in the model that whenever, the main unit fails, it undergoes repair and the 
spare unit starts functioning. After the repair of the main unit, it begins 
operation whereas the spare unit acts as standby. However, if the spare unit 
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fails before the repair of the main unit, it is not repaired and it is scrapped, 
the other spare unit starts operation. In this way, the first failure occurs when 
the main unit is under repair and the last spare unit fails. In the analysis of 
results to obtain total expected cost C(t) in (0,t], they used following costs : 

Cj : cost of repairing the failed unit 

Cj : cost of replacing spare unit. 

Then, C, the expected cost per unit time for an infinite time span is 
C(t) 

C = lim — ^ = Cj. . [expected number of failed units per unit time] 

t -► 00 t 

+ Cj . [expected number of spares per unit time]. 

(6) Okumoto and Osaki (1977) discussed a cost structure in the context 
of obtaining preventive maintenance policies for a standby system. They 
assumed in their model that preventive maintenance is allowed for an operating 
unit only and when it is supported by a standby unit. But the failed unit is 
replaced as it carmot be repaired. 

They have considered the following costs in their model. 

C : replacement cost of a failed unit 

a : preventive maintenance cost per unit time. 

Then, the expected cost C(t) per unit time according to them is given 
by 

C(t) = a . [the limiting probability that a unit is under preventive 
maintenance] + C . [expected number of replacements per unit 
time]. 
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(7) Cleroux and Hanscom (1974) have taken a cost structure which 
takes into account adjustment costs, depreciation costs or interest charges 
which are paid at fixed equidistant intervals of time. They minimized the 
average cost per unit time over an infinite time span to determine optimal 
age replacement policy for a single unit system. 

In notation let 

Cj : cost for each replacement of a failed unit 

C 2 : cost for replacement of a non failed unit 

C 3 (ik) : cost which is incurred at the age ik; i = l, 2.... 

Above costs are known. Further, the sequence of costs [CjCik)] helps 
in keeping the unit operative. The cost [C 3 (ik)] may include adjustment costs, 
depreciation costs, interest charges etc. 

(8) In the problems related to detection of failures as and when they 
occur, some interesting cost structures have appeared in the literature Kumar 
(1975), Gross (1970), Kumar and Kapoor (1981b), Takami et al. (1978)]. 

Kumar (1975) considered the problem of detection of failures in an 
n-component series system. He obtained a condition under which detection 
mechanism is economically feasible using the following cost structure 

Pjj(c) : Prob. [detecting the component as failed whereas infact i*** 
component has failed] 

where c is the amount of money spent on the detection mechanism that 
enhances the probability of correctly detecting the failed component. 


24 


: extra incurred cost if i*** component has failed but failure is 
assigned to component erroneously. 


<^0 


These pij(c) satisfy the following properties. 

n 

2 Pij(c) = 1 for i = 1, 2, ...., n 
j=i 

0 < Pij(c) < 1 for i, j = 1, 2, ..., n 
Pij(O) = ^ for i, j = 1, 2, ..., n 

P;j(c) > ^ for i = j = 1, 2, ..., n and c 5^ 0 

pjj(c) ^ ^ for i j = 1, 2, ..., n and c 0 

Pii(c) t c for i = 1, 2, ..., n 

Pij(c) i c for i j = 1, 2, ..., n 

limpij(c) = dj. 

C-H^OO 

where is Kronecker’s delta 

a.. = 1 1 if ‘ = J 
■J [ 0 if i ^ j 

Gross (1970) is a particular case of the above model when the number 
of components is 2. 

(9) Kumar and Kapoor (1981b) have developed a classification rule for 
a 2-unit series system based on the life-time distribution of the system which 
minimizes the expected cost due to misclassification. 
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They have considered the following costs. 

Cj 2 : cost of classifying the unit 2 as failed whereas unit 1 has failed. 

C 21 : cost of classifying the unit 1 as failed whereas unit 2 has failed. 

(10) Takami et al. (1978) have obtained the solution for the problem 
of optimal allocation of fault detectors in a series system. They have considered 
following costs for their model. 

Cj : cost of the detector allocated to component i 

Cj : loss per unit time caused by the system failure. 

(11) The following reward (cost) structure due to Howard (1971) has 
been extensively applied to obtain expected profit in redundant systems 
operating under different operational conditions, Kumar (1976, 1977c). 

The cost structure is 

Yj : earning rate of the system while it is in the state Sj 

T- : fixed transition cost incurred at the time of transition from 

S; to Sj. 

It may be remarked that the above cost structure allows both types of 
costs viz., fixed and variable costs. 

1.5 SUMMARY OF THE THESIS 

Contribution and Distribution of Chapters TTie review given in the 
preceding pages reveals that for non-maintained systems problem of redundancy 
allocation is important whereas for efficient running of maintained systems it 
is worthwhile to determine optimal failure/repair rates. If a system is in the 
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working state (i.e., it is already designed) one will be interested in having a 
control over the repair process because failure properties might have already 
been prespecified in design stage. So, for maintained systems there is a genuine 
demand of research investigations in the direction of repair/maintenance which 
can help a maintenance engineer/analyst in making an optimal choice of repair 
rule (for a fixed failure rate) and/or in having better understanding of system’s 
economics by studying parameters’ effectiveness from cost considerations. 

Keeping the above in view, the efforts in the thesis are motivated in 
the following directions. 

(i) To provide a decision maker with readily available values of 
repair rates depending upon the earning rates and failure 
characteristics of the system which optimizes steady-state expected 
profit of the system and to provide basic cost equations for the 
better understanding of system’s economic behaviour. 

(ii) To suggest some guidelines to the management regarding 
installation of a new system/or replacing the existing one with 
a new one whose parameters like failure rate, repair rate, earnings 
in various states (viz., operable and failed) are known. 

(iii) Further, in maintained systems, it is of great importance to 
analyze failures i.e., to locate/detect a failure, to classify a failure, 
to isolate the faulty component/components. In this direction, 
two problems of failure detection and classification have also 
been formulated and solved. 

(iv) In the context of the non-maintained systems, a reliability 
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maximization problem in an n-component series system under 
constraint on cost has been discussed. 

While formulating the problems at (i) and (ii) and discussing the 
solution procedures, the cost structure of Howard (1971) has been superimposed 
on the Markov process/semi-Markov process generated by a system model. 
The expected profit thus obtained has been maximized and optimal repair 
rates are determined. This cost structure allows different cost rates (variables) 
in each state and different transition costs incurred at the time when the 
system makes a transition from one state to another. 

For the problem of failure detection in a series system, a cost structure 
involving probabilities of correct detection (depending upon cost spent on 
detection) has been mixed with Howard’s reward structure (1971) and is 
superimposed on a 2-unit series system. Optimal allocation of cost on detection 
mechanisms has been determined when the total cost for detection is fixed. 
Further, a problem of failure classification in a series system has been discussed 
after taking into account costs of misclassification whenever the system fails. 
A decision rule based on the life-time distribution of the system which 
minimizes the expected misclassification cost has been suggested. A procedure 
has also been suggested to examine the discriminatory power of the rule. 

For the problem of cost allocation to different components in a series 
system under a constraint on the total cost, a useful and interesting reliability 
cost function is devised and a solution procedure is developed based on 
dynamic programming. 
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The thesis consists of four chapters : 


Chapter I : 

Chapter 11 : 

Chapter III : 

Chapter IV : 


Introduction 

Optimal Repair Rate in Some Maintained Systems 

Optimal Failure Analysis in Series Systems 

Evaluation of Series Systems Through Dynamic 
Programming Technique. 


Chapter II discusses the profit evaluation studies for some systems with 
repair maintenance. Three models describe single unit maintained systems 
under different operational conditions and one model deals with 2-unit system 
with parallel redundancy. The concept of minor repair has also been applied 
in one of the models. Optimal repair rates maximizing expected profit have 
been obtained analytically for almost all models. The four models on profit 
evaluation studies in this chapter are discussed in the following sections : 


2.2- Profit evaluation in a single unit system 

2.3- Profit evaluation in a single unit system with a less productive 
state 

2.4- Profit considerations in a single unit system with minor repair 

2.5- Profit evaluation in a 2-unit parallel redundant system. 

Chapter III deals with failure analysis in series systems. Problems of 
failure classification and resource allocation on failure detection mechanisms 
have been considered in the following two sections : 

3.2- Optimal failure classification in an n-unit series system. 
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3.3- Optimal allocation of cost to detectors in a 2-unit series system. 

The last chapter discusses application of dynamic programming technique 
to solve problem in series systems. Problem of cost allocation has been 
formulated and solved through dynamic programming technique. This is 
discussed in the following section : 

4.2- Optimum cost allocation in a series system under a constraint 
on the cost. 


30 


CHAPTER II 


OPTIMAL REPAIR RATE IN SOME 


MAINTAINED SYSTEMS 


2.1 INTRODUCTION 


Maintained systems are quite prevalent and useful in modern 
defence/industrial systems. During last three decades a large number of papers 
have appeared to analyse such systems under different conditions. These 
generalizations (contributions) are motivated to the following objectives : 

(i) To generalize the underlying distributions. These distributions 
may be for failure time, repair-time and/or for other features of 
the system units [Nakagawa and Osaki (1974a), Srinivasan and 
Gopalan (1973a), Garg (1963b)] etc. 

(ii) To incorporate certain new parameters in existing models. These 
new parameters may be of interest to system designers and 
maintenance engineers [Kumar and Lai (1979), Arnold (1973), 
Ghow (1975), Gopalan and Saxena (1977), Khalil and Bougas 
(1975), Kapoor and Kapur (1975a), Kumar and Jain (1977)] etc. 

(iii) To obtain and analyse new system parameters which may serve 
useful purpose while making evaluation studies [Kumar (1976), 
Mine and Kawai (1974)] etc. 

Garg (1963b) has taken a general distribution for repair-time and 
analysed a redundant system model. For quite a long time, workers have been 
assuming at least one of the two distributions : failure-time and repair-time 
as exponential. Srinivasan and Gopalan (1973a) have analysed a 2-unit system 
model assuming both the distributions as general. They obtained reliability 
and availability using regenerative properties. Later, Nakagawa and Gsaki 
(1974a) developed a unique modification of MRP and discussed in detail 


reliability characteristics of a 2-unit cold standby system. 

As regards inclusion of new parameters of interest to system 
designers/analysts, some of the interesting parameters include ‘delayed repair’ 
Khalil and Bougas (1975), Kapoor and Kapur (1975a), Gopalan and Saxena 
(1977), Kumar and Jain (1977) : ‘Intermittent repair’ and ‘contact failure’ 
Kumar and Lai (1979), ‘switch failure modes’ Nakagawa (1977a), ‘automatic 
recoveiy’ Arnold (1973) etc. 

Further, Kumar (1976) obtained expected profit for 2-unit systems 
operating under different conditions. This parameter has been suggested as 
the measure of system effectiveness and optimal preventive maintenance 
policies are also discussed. Mine and Kawai (1974) considered a 2-unit parallel 
system with good, degraded and failed states and discussed a preventive 
maintenance policy that maximizes the expected profit rate of the system over 
an infinite time span. Mine and Kawai (1975) discussed optimal replacement 
policy for a unit which assumes any one of several Markov states. The problem 
was formulated as a Markov decision process with a modified Policy 
Improvement Routine and evaluation function is the expected cost per unit 
time over an infinite time span. 

Above review reveals that mostly the contributions are concerned with 
determination of system parameters . and/or inclusion of new features of units 
constituting the system. Nevertheless, efforts have been in the direction of 
helping system analysts and decision makers to take optimal decisions as and 
when they are needed. These papers generally provide answer to the following 
broad question : 


What should be the optimal preventive maintenance policy to maintain 
the system so that the overall system performance (availability, reliability, cost, 
profit etc.) is optimised? 


The purpose of the present chapter is two fold : 


(a) To develop cost equations for a maintaned system and obtain 
thereby expected profit of the system in a given time interval. 

(b) To obtain optimal values for repair rates which may help in 
better management of a maintained system. 


It may be worthwhile to point out here that in all the papers Kumar 
(1976, 1977a, b, c, d) derived expected profit by direct substitution of system 
parameters in Howard’s (1971) expression for g (expected profit). Whereas 
we write basic equations for two models to derive the expression for g. These 
equations help in better understanding of the system’s economic behaviour 
and provide better insight into solution method. However, the result for the 
other two models have been derived by using Howard’s (1971) technique. 

We discuss following models in the subsequent sections : 

Model 1 Profit evalution in a single unit system 

Model 2 Profit evaluation in a single unit system with a less 

productive state 


Model 3 Profit considerations in a single unit system vnth minor 
repair 

Model 4 Profit evaluation in a 2-unit parallel redundant system. 
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22 Model 1 PROFIT EVALUATION IN A SINGLE UNIT SYSTEM 
System Model 

(i) 'niere is a 1-unit system. 

(ii) Failure-time distribution of the system is exponential with 

parameter yS, 

(iii) Whenever the system fails, it undergoes repair. 

(iv) Repair-time distribution of the system is exponential with 

parameter ii. 

(v) Repair facility always exists. 

(vi) After the repair, the system acts as new one. 

(vii) All probability distributions are s-independent. 

(viii) The system earns at a fixed rate so long as it is operable. 

(ix) The system looses at a fixed rate for the time it is under repair. 

(x) There is a fixed transition reward (cost) whenever the system 
changes its state viz., operable (failed) to failed (operable). 

System states and transitions 

The system can be in one of the following 2 states : 

0 : operating 

1 : failed 

To identify the system at any instant, let 
Sq : system is operating 
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Sj : system is in failed state and is under repair. 

Notation 

/S : constant failure rate of the system 

^ : constant repair rate of the system 

rjj : cost per unit time or earning rate in (i = 0, 1) 

ry : transition cost for a transition from state Sj to Sj, 

(t, j = 0, 1) 

Vi(t) : total expected earning of the system for a future period of 
length ‘t’ given that at t = 0 the system was in Sj (i = 0, 1) 

A : small interval 

f(s) : Laplace transform of f(t), i.e., 

f(s) = / e-*‘ f(t) dt, Re (s) > 0. 

0 

Cost Equations and System Profit 

Following Howard (1960), the following cost equations can be written 

Vo(t+A) = [rooA+Vo(t)][l-M]+[roi+Vi(t)]^A (2.2.1) 

V,(t+A) = [r,jA+V,(t)][l-A4A]+[r,o+Vo(t)luA (2.2.2) 

Initially V; (0) = 0 (i = 0, 1). 

The equation 2.2.1 is interpreted in the following manner. 

During the short time interval A, the system which started in ‘0’ state 
at t = 0 may remain in ‘0’ state or make a transition to state T’. If it remains 


in ‘0’ state for time A, the system’s earning will be rpo A plus the expected 
earning it will make in the remaining time of ‘t’ units viz., Vj,(t). The probability 
that it remains in state ‘0’ for time A is 1 minus the probability that it makes 
a transition to state T’ in A, i.e., 1-/5A . In case the system makes a transition 
to state T’ during time interval A, the system’s earning will be r^j plus the 
expected earning it will make when the system starts in state T’ with time t 
remaining viz., Vj(t). The sum of the products of earnings with their respective 
probabilities gives the total expected earning. Equation 2.2.2 can be interpreted 
in the same manner. 

Expanding equations 2.2.1 and 2.2.2 we get 


Vo(t+A) = rooA+Vo(t)-^Vo(t)A+/3rojA+^Vj(t)A+0(A) (2.2.3) 

Vi(t+A) = riiA+Vi(t)-/zVi(t)A+/^rioA+/iVo(t)A+0(A) (2.2.4) 

Neglecting terms of higher order of A and letting A 0, we obtain 
the following system of differential equations 




Vi(t) = rii+/^rio+/^Vo(t) 


(2.2.5) 

( 2 . 2 . 6 ) 


Applying Laplace transforms to the set of equations 2.2.5 and 2.2.6 and 
employing initial conditions Vi(0) = 0, (i = 0, 1) we get 


T7 , . Too+^^oi . iS ^ 


(2.2.7) 


s(s+/^)^s+/^^o^^^ 


( 2 . 2 . 8 ) 


On substituting the value of Vj(s) form (2.2.8) in (2.2.7), we get 


Vo(s) 


(roo+^roi)s+i“(roo+^roi)+/^(rii+/^rio) 

s^(s+yS+/i) 


(2.2.9) 


On substituting the value of Vo(s) from (2.2.7) in (2.2.8), we get 


Vj(s) = 


s\s+^+fi) 


( 2 . 2 . 10 ) 


On inverting (2.2.9) and (2.2.10), we get 


-yS(r„+/ir,o-roo-^ro,)e-(^^^>‘] (2.2.11) 


^l(0 ~ j3+/W 

- /.(roo+Aoi-rii-/^r,o)e-'^"^)‘] (2.2.12) 


In order to exhibit transient behaviour of the system and thus usage 
of results arrived above in (2.2.11) and (2.2.12), we give below an illustration. 


Illustration 

Let us specify various parameters for the purpose of numerical 
illustration. We give the results in the subsequent tables (2.2.1-2.2.4) for 
different values of y3, pi, r^o, r^, r^j and rjg 
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Define VjCt) = V„(t + 1) - V„{t) 
and V5(t) = V,{t + I) - Vj(t) 

From table 2.2.1, it may be observed if the system starts from operable 
state at t = 0 it achieves steady-state after 13 units of time and the steady-state 
expected profit is 3.5238. If the system starts fi’om failed state it achieves 
steady-state after 14 units of time and the steady-state expected profit is 
3.5238. 


t 

Vo(t) 

vi(t) 

Vi(t) 


0 

0 

13.2600 

0 

-27.6320 

1 

13.2600 

6.9308 

-27.6320 

-7.3788 

2 

20.1908 

4.7161 

-35.0108 

-0.2914 

3 

24.9069 

3.9410 

-35.3022 

2.1887 

4 

28.8479 

3.6698 

-33.1135 

3.0567 

5 

32.5177 

3.5749 

-30.0568 

3.3603 

6 

36.0926 

3.5417 

-26.6965 

3.4666 

7 

39.6343 

3.5301 

-23.2299 

3.5038 

8 

43.1644 

3.5260 

-19.7261 

3.5168 

9 

46.6904 

3.5246 

-16.2093 

3.5213 

10 

50.2150 

3.5240 

-12.6880 

3.5230 

11 

53.7390 

3.5239 

-9.1650 

3.5235 

12 

57.2629 

3.5239 

-5.6415 

3.5237 

13 

60.7868 

3.5238 

-2.1178 

3.5237 

14 

64.3106 

3.5238 

1.4059 

3.5238 

15 

67.8344 

3.5238 

4.9297 

3.5238 

16 

71.3582 

3.5238 

8.4535 

3.5238 

17 

74.8820 


11.9773 



/S = 0.25, jU= 0.8, Too = 20.0, rji = -50.0, roi = -3.0, rjo=4.0 


Table 2,2.1 


It may be observed that initially there is steep fall in Vo(t), whereas, 
as the time progresses, this decrease reduces and stabilizes to a fixed value 
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which is called steady-state ejqpected profit of the system, Kumar (1976), 
Further, in the beginning Vj(t) increases rapidly and then it stabilizes to the 
same fixed value viz., expected steady-state profit. 

From table 2.2.2, it may be observed that steady-state is achieved after 
11 units of time if system starts in operable state and the expected profit is 
5.9048. Further, it may be noted that when earning in the failed state increases 
(cost spent on repairs decreases), expected profit increases and a reduction 


of 20% in repairs ensures 67.5% increase in the expected profit. 


t 

Vo(t) 

V|)(t) 

Vi(t) 

v!,(t) 

0 

0 

14.1668 

0 

-20.5340 

1 

14.1668 

8.7960 

-20.5340 

-3.3472 

2 

22.9628 

6.9166 

-23.8812 

2.6672 

3 

29.8794 

6.1982 

-21.2140 

4.7718 

4 

36.1382 

6.0286 

-16.4422 

5.5083 

5 

42.1668 

5.9481 

-10.9339 

5.7660 

6 

48.1149 

5.9200 

-5.1679 

5.8561 

7 

54.0349 

5.9100 

0.6882 

5.8878 

8 

59.9449 

5.9067 

6.5760 

5.8988 

9 

65.8516 

5.9054 

12.4748 

5.9027 

10 

71.7570 

5.9050 

18.3775 

5.9040 

11 

77.6620 

5.9048 

24.2815 

5.9045 

12 

83.5668 

5.9048 

30.1860 

5.9047 

13 

89.4716 

5.9048 

36.0907 

5.9047 

14 

95.3764 

5.9047 

41.9954 

5.9048 

15 

101.2811 

5.9048 

47.9002 

5.9047 

16 

107.1859 

5.9048 

53.8049 

5.9048 

17 

113.0907 

5.9047 

59.7097 

5.9048 

18 

118.9954 


65.6145 



= 0.25, 0.8, roo=20.0, rjj = -40.0, rQj=-3.0, rj(, = 4.0 


Table 2.2 JS 
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! 

Vo(t) 

V!,(t) 

v,(t) 

v'lO) 

0 

0 

17.8065 

0 

-26.1809 

1 

17.8065 

10.9983 

-26.1809 

-4.3946 

2 

28.8048 

8.6159 

-30.5755 

3.2293 

3 

37.4207 

7.7821 

-27.3462 

5.8971 

4 

45.2028 

7.4904 


6.8308 

5 

52.6932 

7.3883 

-14.6183 


6 

60.0815 

7.3525 

-7.4608 


7 

67.4340 

7.3401 

-0.1890 

7.3117 

8 

74.7741 

7.3357 

7.1227 

7.3258 

9 

82.1098 

7.3342 

14.4485 

7.3307 

10 

89.4440 

7.3336 

21.7792 

7.3324 

11 

96.7776 

7.3334 

29.1116 

7.3330 

12 

104.1110 

7.3334 

36.4446 

7,3332 

13 

111.4444 

7.3333 

43.7778 

7.3333 

14 

118.7777 

7.3334 

51.1111 

7.3333 

15 

126.1111 

7.3333 

58.4444 

7.3333 

16 

133.4444 

7.3333 

65.7777 

7.3334 

I 17 

140.7777 


73.1111 



P = 0.25, fJL = 0.8, Too = 25.0, rjj = -50.0, Fqj = -3.0, Tjo = 4.0 


Table 22 J 

In this case, steady-state is achieved at t = 13. Further, it may be noted 
that an increase of 25% in the earning rate of the system while it is operative, 
there is 108.1% increase in steady-state profit of the system. 
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0 

0 

13.6488 

0 

- 23.5955 

1 

13.6488 

8.3342 

- 23.5955 

- 2.3366 

2 

21.9830 

6.8114 

- 25.9321 

3.7542 

3 

28.7944 

6.3752 

- 22.1779 

5.4991 

4 

35.1696 

6.2502 

- 16.6788 

5.9995 

5 

41.4198 

6.2144 

- 10.6793 

6.1424 

6 : 

47.6342 

6.2041 

- 4.5369 

6.1835 

7 

53.8383 

6.2012 

1.6466 

6.1953 

8 

60.0395 

6.2003 

7.8419 

6.1986 

9 

66.2398 

6.2001 

14.0405 

6.1996 

10 

72.4399 

6.2000 

20.2401 

6.1999 

11 

78.6399 

6.2001 

26.4400 

6.2000 

12 

84.8400 

6.2000 

32.6400 

6.2000 

13 

91.0400 

6.2000 

38.8400 

6.2000 

14 

97.2400 

6.2000 

45.0400 

6.2000 

15 

103.4400 


51.2400 



yS = 0.25, /^ = 1.0, rQQ = 20.0, rij = -50.0, rQj=-3.0, rjg = 4.0 


Table 22.4 

Table 2.2.4 shows that faster repair results in increase in profit. An 
increase in repair rate by 25% gives an increase of about 76% in steady-state 
profit of the system. 

We now depict transient behaviour of the economics of the system so 
as to bring into light some interesting facts which may be helpful for systems 
which are operated for a short time span. Behaviour of Vo(t) with respect to 
variation in t has been given in Fig. 2.2.1. Some of the observations are as 
under : 

(i) The total earning of the system increases initially. As time 
progresses, the earning reaches a maximum limit and then starts 
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FIGURE; 2-21 

GRAPH SHOWING THE BEHAVIOUR OF TOTAL 
EXPECTED EARNING Vo (t) AGAINST ‘t* 

(WHEN tIS NOT LARGE ENOUGH) 
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PROFIT REGION 


7 8 9 

LOSS REGION 




decreasing. At first sight it may look ridiculous that earning of 
the system in time t decreases as the time interval t increases. 
But, a close look in the system structure parameters viz., failure 
rate, repair rate etc. reveals that the earning does not solely 
depend upon the time but it is influenced by several other factors 
including failure and repair rates of the system. When the earning 
in a larger interval reduces, this merely implies that a major 
share of the money is spent on repairs of the system. 

(ii) Further, just after 6 units of time, the system starts running in 
loss and it continues to be in this region in steady-state. More 
elaborately, after initial adjustments, the system runs in loss. This 
happens because of the choice of cost parameters and system 

parameters. 

Expected Profit Optimization 

The expected profit per unit time denoted by g from (2.2.11) or (2.2.12) 
after allowing t tend to infinity can be easily obtained. Thus, we get 

(rm+^r„,)u+(rii+/^rio)^ ('2.2.131 

g = p+Ji 

The expected profit g, thus depends upon failure rate repair rate 
fi and r-j (i, j = 0, 1). It may be observed that it will be of interest to obtain 
the value for the repair rate which maximizes the expected profit when r^ 
and ^ are given (fixed). For the purpose, let us consider 

= = = (2.2.14) 
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Further, if we assume 


Cj > 0, ^ > 0 

A<0 
B>0 

rji = f 4 (/^) = a+c^+c^\ C 2 ,C 3 < 0 , a is a constant 

On making use of (2.2.15) in (2.2.13) the expected profit g takes the 


1-00 = 

^10 ~ ^ 3 ^) ~ 


(2.2.15) 


form 


® fi+d 


(2.2.16) 


where 


a=^C3, b=j + ^(A + B + C2), c=a0, d=^. 


The value of ft which maximizes the expected profit m (2.2.16) is 
obtained by differentiating g in (2.2.16) with respect to p. and equating it to 


zero i.e., 


dg _ a|U^+2ad^+bd-c _ ^ 
dp {p+df 


^ 2 p^ + 22 idp+hd-c = 0 (2.2.17) 

The optimum value of repair rate p denoted by p* which maximizes 
g denoted by g* is given by 


-?.ad± Vr2ady-4a(bd-cT ^ -ad±V^ 
2a a 


(2.2.18) 














if"!-; 

H 


Wi 
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where 


£y = (ad)^ - a(bd - c) 

Since fi* is, finite and cannot be negative, being repair rate, to decide 


fi* gives maximum we have to further examine the sign of 2 




\i fi* > 0 and ^ ^ maximum g i.e., 

the sufficient condition for/f ’>0 to be absolute maximum point 
for g is 

+ ^^(A + B + C2) + - Cj < 0. 


(ii) However, ii pt >0 and - 2 I „ — „* ^ ^ minimum 

v '' d|« I - /♦ 

point for g i.e., when 

+ ^^(A+B+C2) + ajS - Cj > 0 
and absolute maximum is achieved at =0 for which g =«• 

We now proceed to illustrate the use of results in (2.2.18). We present 
optimum repair rate fi against various values of failure rates, ^ for different 
sets of constants a, c^, C 2 , C 3 , A and B (where a = 0) in the Figure 2.2.2. 
Some of observations are : 

(i) It may be observed that in general {i decreases as ^ increases. 
Further the decrease in the value of for small (initial) values 
of p is quite substantial as compared to the rate of decrease in 
the value of for large values of jS. To elaborate this point 


45 



figure: 2-2-2: GRAPH SHOWING THE VARIATION OF OPTIMUM 
REPAIR RATE(^?AGAINST FAILURE RATE C-P) 
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further, we present the results in table 2.2.4a corresponding to 
curve 1 in Figure 2.2.2. From the same table we observe that 
the expected profit g* also behaves in similar manner. Behaviour 


of other curves also follow the same pattern. 



« 

lOOSISEBI 

★ 

g 

g’03)-g*03 + .1) 

0.05 

6.27 

2.80 

196.47 

136.21 

0.15 

3.47 

0.95 

60.26 

28.56 

0.25 

2.52 

0.57 

31.70 

12.74 

0.35 

1.95 

0.42 

18.96 

7.23 

0.45 

1.53 

0.33 

11.73 

4.54 

0.55 

1.20 

0.28 

7.19 

2.97 

0.65 

0.92 

0.25 

4.22 

1.95 

0.75 

0.67 

0.21 

2.27 

1.21 

0.85 

0.46 

0.19 

1.06 

0.70 

0.95 

0.27 

0.17 

0.36 

0.31 

1.05 

0.10 


0.05 



Table 22.4a 


(ii) With the increase in Cj, the value of fi increases but shows a 


steep increase for smaller values of for example 



fj* ixom curve (1) 

/<*from curve (2) 

Difference 

0.25 

2.51 

3.16 

.0.65 

1.00 

0.18 

0.55 

0.37 


(iii) With the increase in Cj, the value of fi first increases slowly 


and then significantly, for example 



IIIIIQI^^ 


Difference 

0.25 

2.51 

2.55 

0.04 

1.00 

0.18 

0.48 

0.30 


(iv) With the increase in Cj, the value of fi first increases significantly 
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and then slowly, for example 



/i^from curve (1) 

^*from curve (4) 

Difference 

0.25 

2.51 

3.30 

0.79 

1.00 

0.18 

0.29 

0.11 

(v) With 

the increase in A, the value of first increases slowly 

and then significantly, for example 




fi* from curve (6) 

Difference 

0.25 

2.51 


0.03 

1.00 

0.18 


0.16 


Further, if one is confronted with the problem of making a choice 
regarding installing of a production system whose various parameters like 
failure rate yS, repair rate n and the earning rates in operable state and failed 
state, transition costs for making a transition from operable state to failed 
state and vice versa i.e., r^ are specified. Then one should choose a system 
which gives more expected profit than others making use of (2.2.13). In what 
follows we present the table 2.2.5 for given system under consideration. 


System 

Parameter 

1 

2 

3 

4 

5 

6 

7 

8 

§ 

0.5 

0.6 

0.7 

1.0 

0.5 

0.8 


1.0 


1.5 

1.3 

1.8 

1.5 

1.5 

2.0 


2.5 

*^00 

20.0 

25.0 

30.0 

40.0 

35.0 

50.0 

55.0 

55.0 

^'ll 

-30.0 

-30.0 

-35.0 

-50.0 

-40.0 

-65.0 

-70.0 

-75.0 

^01 

-4.0 

-4.0 

-2.0 

-2.0 

-3.0 

-2.0 

-3.0 

-4.0 

*^10 

3.0 

3.0 

3.0 

3.0 

4.0 

5.0 

4.0 

2.0 


7.12 

^ 7.22 

12.30 

4.60 

16.62 

1 18.85 

14.00 

16.42 


Table 2.2.5 

Thus we find from the above table that the expected profit g is maximum 
for the 6*** system and therefore, management should opt for this system. 
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23 Model 2 : PROFIT EVALUATION IN A SINGLE UNIT SYSTEM 
WITH A LESS PRODUCTIVE STATE 

In this section we discuss a model of a single unit consisting of two 
modules. The failure of one module brings the system in a less productive 
state whereas the failure in the other module brings the system to failed state. 

System Model 

(i) There is 1-unit system consisting of type I module and type 11 
module. 

(ii) From operable state, system goes to failed state on the failure 
of type I module. 

(iii) From operable state, system goes to less productive state on the 
failure of type II module. 

(iv) From less productive state system goes to failed state on the 
failure of type I module. 

(v) In the failed state, repairs are carried out to bring the system 
back to operable state. 

(vi) Failure time distributions of type I and type II modules are 
assumed to be exponential. 

(vii) Repair time distribution is assumed to be general. 

(viii) In failed state, if there is any failed module of type II, it is 
repaired free of cost together with repair of type I module. 

(ix) After repair, system acts as new one. 
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(x) All probability distributions are assumed to be s-independent. 

(xi) The system earns (loses) at a fixed rate in each state which can 
be different for each state. There is a fixed transition reward 
(cost) whenever the system changes its state. 

Define the following system states to identify the system at any instant. 

System states and Transitions 

Sq : Operable or fully productive, 

Sj : Less productive, type II module is in failed state. 

S 2 : Failed state, type I module or type I and type 11 modules 

are in failed state and the system is under repair. 

The system is up in Sq and Sj, and it is down in S 2 . The transitions 
between states are given in Figure 2.3.1. 



Figure 2.3.1. Transitions between various states 
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NOTATION 

/Sj : constant failure rate of Type I module. 

^2 ■ constant failure rate of Type II module. 

g(t) : pdf of repair time 
fi : expected time for repair 

Tj : mean unconditional waiting time in Sj (i = 0, 1, 2) 

pjj : one step transitional probability from Sj to Sj 

jr, : steady-state probability in Sj (i = 0, 1, 2) 

rjj : transition cost for a transition from Sj to Sj (for i ^ j). 

rjj : cost per unit time or earning rate in Sj 

g : total expected earning per unit time in steady-state or 

steady-state expected profit. 


Expected Profit of the System 

Observing state transitions from figure 2.3.1, various state transition 
probabilities can be written as : 


Poi = 

Po2 = dt = 

00 

P20 = J' g(t)dt=l 

0 

Pi2 = /A dt=l 


(2.3.1) 
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Further, it is easy to see that 



e“^2‘ dt 


KW2 


Tj = /e"'^>‘ dt = i 
0 1^1 



(2.3.2) 


Steady-state probabilities n- can be obtained by following relations of 
Howard (1971) and from (2.3.1) 

JTj = 2 Pi/Tj (2.3.3) 

i 

and 

2 Jt, = 1 (2.3.4) 

i 

Therefore, 


_ , _ ^2 , _ 


(2.3.5) 


Following Howard (1971) expected profit g can be written as 


2 JiiTa 


(2.3.6) 


where 


qi = rii + ST 2 

j 


m 
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Substituting values of from (2.3.1), Tj from (2.3.2) and jT; (2.3.5) in 
(2.3.6) and after simplification g can be obtained as 


(^00 "^^1^02) 1 1 '*’^1*^12) 1 1 C^22 ^20) 


(2.3.7) 


Profit Optimization 

Let the repair rate n of the system is controllable variable, expected 
profit, g can be optimized for two cases as under : 


Case I 

When the earning in the failed state is a continuous function of repair 
rate. Let us specify tjj and rjj as follows : 


roo = A, A > 0, 


Tjj = B, A > B > 0 

% = , a<0; Cj, C 2 < 0 


and 

r-; = 0 for i 9 ^ j 

Ij J 


(2.3.8) 


Substituting these value of Tqo, r^, r 22 from (2.3.8) in (2.3.7) and after 
simplification we get 


_ K^^+I^+N 
^ P/i+Y 


(2.3.9) 


where 
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K=(^,+^^A 

L=^,A+^JB+C8,+/3j)e,C, 

P = (fii+Pi) 

Y = 


(2.3.10) 


To determine the value of n which maximizes g for given 
known values of cost parameters, we differentiate g in (2.3.9) with respect to 

H and equate it to zero i.e., ^ = 0, which gives the value of fi denoted by 


maximizing g denoted by g*, so 




KY±VkV-KP(LY-NP) 

KP 


(2.3.11) 


Since n' is finite and cannot be negative, being repair rate, to decide 


* d^g 

fi gives maximum, we have to further examine the sign of 


(i) 


(ii) 


fi =fi 


If /u* > 0 and 


fi =fi 


, < 0, fi' provides absolute maximum 


g* i.e., the sufficient condition for jU > 0 to be absolute maximum 
point for g is < A^i+E ^2 


However, if /f ’ > 0 and 


fi =fi 


* > Q, fi is the minimum 


point i.e., when 03j+^2)(/^iQj~ACj+a) > 2 ^j+E ^2 . is the 
minimum point and absolute maximum is achieved st ft = 0 
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for which g* = a. 


Illustration 

In what follows we shall illustrate the results with the help of numerical 
example. Let us consider a system for which the various values of parameters 
are given as below : 

= A = 50.0, a = -50.0, € 2 = -2.0 

y32 = 0.2, B = 40.0, Ci = ^.0 

On substituting these values A, B, a, Cj and C 2 (2.3.10) and 

(2.3.11), we obtain /f’ = 6.916 

For /4’ = 6.916, g’ can be obtained from (2.3.9) i.e., 

g* = 40.245. 

Case II 

When the earning in the failed state is a discrete function of repair 

rate. 

In many practical situations the earning of the system in failed state 
i.e., r 22 is not a continuous function of the repair rate fi and in fact r 22 is 
discrete function of repair rate fi. 

Let the management have the option to choose a repair policy out of 
a number of given repair policies (repair rates) and their corresponding earning 
rates (r 22 ) in failed states. 

In such cases, g is computed for each fi and corresponding r 22 from 
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(2.3.7) and that repair policy is chosen which gives more expected profit than 
others. 

In case more then one repair policy gives the same amount of expected 
profit, any one of them can be selected. 

Illustration 

In what follows we shall illustrate with a numerical example. Let the 
values of the various parameters of the system under consideration are specified 
as below : 

Too - 50,0, )8j = 0. 1 , r^j = -2.0, r jj = -2,0 
rjj = 40 . 0 , ;32 = 0 . 2 , ro2 = - 3 . 0 , r2o = 4 . 0 . 

Then the values of g computed from (2.3.7) for known values of repair 


rates, and r 22 are presented in the table 2.3.1. 



1.0 

2.0 

3.0 

1 

4.0 


5.0 

7.0 

r 22 

- 60.0 

- 68.0 



- 88.0 

- 105.0 

- 160.0 

g 

33.969 

38.063 

39.548 

40.235 

40.510 

40.457 

40.502 


Table 23.1 

From the above table we observe that for repair rate, /<=4.5, g is 
maximum and hence the management should follow this repair policy. 


2.4 Model. 3 : PROFIT CONSIDERATIONS IN A SINGLE UNIT 

SYSTEM WITH MINOR REPAIR 

System Model 

(i) There is a l-unit system. The system can go to either failed 
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state or less productive state. 

(ii) When the system goes to less productive state, the management 
can opt : 

(a) To carry out minor repair in order to bring the system back 
to fully productive state. 

(b) To allow the system to fail. 

(iii) In the failed state repairs are carried out to bring the system 
back to operable state. 

(iv) Failure-time distributions namely going from operable to failed 
state and less productive state are assumed to be exponential. 

(v) Repair-time distributions are assumed to be general. 

(vi) All probability distributions are assumed to be s-independent. 

(vii) The system earns (loses) at a iSxed rate in each state which can 
be different for each state. There is a fixed transition reward 
(cost) whenever the system changes its state. 

Define the following system states to identify the system at any instant. 

System states and Transitions 

Sq : Operable or fully productive 

Sj : Less productive 

$2 : Failed and is under repair 

We consider two cases - one under which the system is sent for minor 
repairs as soon as it enters Sj and second when the system is allowed to 


57 



enter 83 in due course of time. This is depicted in transitional diagrams 
(Figure 2.4.1 and 2.4.2). 



NOTATION 

: Constant rate of going from to Sj 

^2 '■ Constant rate of going from Sq to 83 

^3 : Constant rate of going from 8j to 83 

gj(t) : pdf of minor repair time 

g3(t) : pdf of repair time 

: expected time for minor repair 
/ij : expected time for repair 

Tj ; mean unconditional waiting time in 8j (i = 0, 1, 2) 

P;j : one step transitional probability from 8; to Sj 
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jr. : steady-state probability in Sj (i=0, 1, 2) 

Fjj : transition reward for a transition from Sj to Sj 

r^j : earning rate in S; (i = 0, 2) 

: earning rate in Sj when minor repairs are being carried out 

: earning rate in Sj when minor repairs are not being carried 
out 

gj^ : expected profit per unit time when the system is under minor 
repair 

gj^r : expected profit per unit time when minor repairs are not 
carried out in Sj. 


Expected Profit of the System 
Case I 

Observing state transitions from figure 2.4.1, various transition 
probabilities can be written as 


Poi = { 

P02~^ ^2® ® 

00 

Pio=-^ gi(Odt=l 

0 

00 

P20 = J’ g2(t)dt=l 

0 

Further, it is easy to see that 


(2.4.1) 
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(2.4.2) 


00 



0 




f*2 


Also steady-state probabilities can be obtained by following relations 
of Howard (1971) and p^j from (2.4.1) 


= 2 Pij^i 

i 



(2.4.3) 

2 = 1 

i 



(2.4.4) 

Therefore, 




1 



(2.4.5) 

^0 ~ 2 ’ 

2(^,+^,) ' 

^ 2(P^+^) 


Following Howard (1971), expected profit per unit time, g can be 
written as 


2 Tj qi 

g = V tt" 

2 Tj 

i 

where qi = ry + ^E Py ^ 
j 


Substituting values of Pjj from (2.4.1), T; from (2.4.2) and from (2.4.5) 
in (2.4.6) and after simplification, gj^, can be obtained as 
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_ ^2(^11 ~^/^l^lo)'^^2^l(^ 22~*~'^2^2o) 

1 ^ 2 '*’^ 1^2 "^^ 2^1 

This can also be written as 


/2 jA+B 
"" fip+E 


(2.4.7) 


where A =/22(roo+^iroi+te2)+^i^2rio+M2+^2/^2r2o 

D = ^2 ■^-“2 
E = 



Case II 

Observing state transitions from figure 2.4.2, various transition 
probabilities can be written as 
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(2.4.8) 


0 

00 

P20= / g2(0<it=l 
0 

Further, it is easy to see that 
To = / dt=ra 


1 ' 1'2 


Tj = /e'^3‘ dt=^ 
0 Fa 


T 


(2.4.9) 


Also steady-state probabilities jTj can be obtained by making use of 
(2.4.3), (2.4.4) and p^- from (2.4.8). 


Jtfi 


A+i ^2 


jr, = 




JTn ~ 


A +^2 


3 ^ 9 , + 2^2 ’ 2 


(2.4.10) 


Substituting these values of p^ from (2.4.8), Tj from (2.4.9) and jTj from 
(2.4.10) in (2.4.6) and after simplification we get, 


Sn 


/i^3-l-y3ja2+A0?i +y^2) 


This can also be written as 


(2.4.11) 



where 


G = P-J^2^n 
H = ^i+^2+/“2 
E = 

Decision Criterion Regarding Profit 

Now it is evident that when > g^., the management of the system 
should follow procedures to carry out minor repairs for acquiring more profit. 

When gj^j = gpj, it is immaterial whether minor repairs are carried out 
in less productive state or the system is allowed to fail. 

When gj^ < gj^, it is advisable, not to follow minor repair procedures. 

For gj^ > gj., 

Then, from (2.4.7) and (2.4.11) 

/fjA+B ^jF+G 
yUjD+E ~ ySjH+E 

=> /fj/S3(AH-FD)+/^j(EA-GD)+^3(BH-FE)>GE-BE (2.4.12) 

Let 



Therefore, (2.4.12) reduce to 

K[/S^(AH-FD)+^3(EA-GD)]>GE-BE-^3(BH-FE). 
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Then, we get critical value of K denoted by K‘ as 
(GE-BE)-A(BH-FE) 

K = -7 — (2.4.13) 

yS^(AH-FD)+)$3(EA-GD) ^ ^ 

For this value of K’, = gj^,. 

Further, notice if K < K*, it is profitable not to carry out minor 
repairs, since this leads to less profit. 

However, when K > K*, it is profitable to carry out minor repairs. 


ILLUSTRATION 

Let us consider a system for which the various values of parameters 


are given as below 

except /Uj and ySg 


A =0.6, 

Too = 10.0, 

r^ = -5.0, 

ri2 = -3.0 

A = 0.4, 

^01 ~ '2.0, 

rrj = 5.0, 

0 

in 

II 

0 

/'2=3.0, 

o' 

T 

II 

rjo=3.0, 

r22 = -15.0. 


Substituting these values in (2.4.13), the critical values of K* are given 
in Table 2.4.1 for various values of 


^3 


k’ 

0.1 

22.505 

0.2 

10.058 

0.3 

6.211 

0.4 

4.406 

0.5 

3.379 

1.0 

1.502 

1.5 

0.949 

2.0 

0.690 

2.5 

0.541 


Table 2.4.1 
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Now suppose in the given problem ^3 = 1.0 and = then, K = 1.0 
which is less than K (K =1.502 from the above table) and therefore, it is 
profitable if minor repairs are not carried out. 

Now suppose /S 3 = 1.0 and ;/j=2.0 then K=2.0 which is greater than 
K' and therefore, it is profitable to carry out minor repair. 

2.5 Model 4 : PROFIT EVALUATION IN A 2-UNIT PARALLEL 
REDUNDANT SYSTEM 

System Model 

(i) The system consists of two units, I and II. Both the units operate 
simultaneously. 

(ii) Failure-time distributions of units I and II are exponential with 
rates /Sj and respectively. 

(iii) Whenever one of the two units fails, the operative unit continues 
to operate whereas the failed unit undergoes repair. 

(iv) On failure of both units, they undergo repair. 

(v) Repair-time distributions of units I and II are exponential with 
rates /fj and respectively. 

(vi) All probability distributions are assumed to be statistically 
independent. 

(vii) After repair, units act as new ones. 

(viii) The system earns (loses) at a fixed rate in each state which can 
be different for each state. There is a fixed transition reward 
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(cost) whenever the system changes its state. 

Define the following system states to identify the system at any instant. 

System States and Transitions 

Sg : Both the units are operating. 

Sj : Unit II is operating while unit I fails and undergoes repair. 

$2 : Unit I is operating while unit II fails and undergoes repair. 

S 3 : Both the units are in failed state and are under repair. 

The system is up in Sg, Sj and S 2 and it is down in S 3 . Transitions 
between various states are given in Figure 2.5.1. 



Figure 2.5.1 : Transitions between various states 
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NOTATION 


: constant failure rate of Unit I. 

^2 • constant failure rate of Unit II. 

//j : constant repair rate of Unit I. 

fi 2 : constant repair rate of Unit II. 

r^j : cost per unit time or earning rate in state Sj. (i = 0, 1, 2, 3). 

rji : transition cost for a transition from state S: to S:, i 5 ^ j 

(i, j = 0, 1, 2, 3). 

Vj(t) : total expected earning of the system for a future period of 
length ‘t’ given that at t = 0, the system was in Sj 
(i = 0, 1, 2, 3). 

Cost Equations and System Profit 

Following Howard (1960), we can write the following cost equations : 


V„(t+A) = [r;^+V,(t)][l-/3,A][l-ftA]+tr„,+V,(t)] 

X [l-ftA13,A+[r<„+V,(t)lII-^,Al^,A 

(2.5.1) 

V,(t+A) = [r„A+V,(t)][l-/^,AJ[l-ftAl+[r,„+V„(t)] 

X [l-/9,Alu,A+[r„+V3(t)][l-/i,A13,A 

(2.5.2) 

VjCt+A) = [r33A+V3(t)][l-//,A)[l-^,A]+[r3;,+V„(t)] 

X [l-ftAlu3A+[r33+V,(t)][l-/<,A)S,A 

(2.5.3) 

V3(t+A) = [r33A+V3(t)][l-//,A](l-ftAJ+[r3,+V,{t)l 

X [l-/<,AlU3A+[r33+V2(t)][l-/<3A]«,A 

(2.5.4) 
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where A is a small interval. 


Initially Vj(0) = 0 for i = 0, 1, 2, 3. (2.5.5) 

Neglecting terms of order of 2 and higher of A, and letting A -»■ 0, we 
obtain following system of differential equations from (2.5.1) to (2.5.4). 

+ &3,+/82)1Vo«) = roo+/S.ro.+to2+AV,(t)+y5jVj{t) (2.5.6) 

I* + (f'i+/52)]V,(t) = r„+//,r,„+/?,r„+/2,V„(t)+^2V3(t) (2.5.7) 

~ ^*22 ’^(^2^20 "^^1^23 ^ 3(0 (2.5.8) 

['^ C^i'^A*2)3^3(^) ~ ^33'*'^2^31'*’'^1^32'^^2^i(0'1’/^1^2(0 (2.5.9) 

Let f (s) denote the Laplace transform of f(t), e.g., 

f(s) = / e"‘'‘f(t)dt, Re (s)>0. 

0 

Taking Laplace transform of (2.5.6) to (2.5.9) and employing initial 
conditions (2.5.5), we get after simplification : 


Vo(s) = 


+ 

ftV,(s)+ftV,(s) 

(2.5.10) 

S(S+^j+i?2) 

(s+^1+^2) 

V:(s) = 


t 

/'.V„(s)+ftV 3 (s) 

(2.5.11) 

s(s+/^l+^2) 

T- 

(s+Ai+^2) 

v^Cs) = 


+ 

/' 2 V„(s)+ftV 3 {s) 

(2.5.12) 

s(s+i?i+/^2) 

(s+^l+/<2) 

¥3(5) = 

r 33 +/' 2 r 3 i+/'ir 32 

1 

2 ^ j(s) +/^ j V 2 ( 5 ) 

(2.5.13) 

S(S+/^l+/^2) 
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Further simplification of (2.5.10) to (2.5.13) gives 


Vo(s) 





(2.5.14) 


where 


Rq ~ ^00 ^1^01 2^02 

Rj =rjj + ■•■^2% 

R2 “ ^22 -“2% ^ 1^23 

R3~^33 ^2% '*■'“1^32 


= (S+yUj+^2)(s+^l+/^2)(S+/«l+/^2)~^]^](S+/“l+^2)“^2/“2(S+/^l+<“2) 
A 2 = (S+^i+/^2)(S+/«1+/^2)-1^i/^1+^2/^2 
A 3 = (S+yS2+/<j)(s+//j+/22)+^j^j-^2/^2 


A 4 = 2 s+^ 1 +/^j+^ 2 +/^ 2 - 
On inverting (2.5.14) we get 


a^e a,e 

Vo(t) = Ro[-^ + 


^ ] + B,[-^ + ^ + 

^4 -^2^5 ajas 84 


+ B2[- 


l-e'V _ e"”!* 


+ 


e"^2‘ 

+ ] 




+ B3[' 


841-83 82^6 ® 2 *-aj^e 


31- „2 


+ 


(2.5.15) 
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where 


~ 

®2 ~ [%0^l'^/*2)’^^1^2J^O*^(^l'*’'^2)^l^l'^(^'^^3^2^2'^^1^2^3 
® 3 ~ 


and 


aj— ^j+/ij, Z^—^2^P'2'> —^1 


Similarly expressions for Vj(t), ¥ 3 ( 1 ) and ¥ 3 ( 1 ) can be obtained. The 
total expected earning of the system per unit time in steady-state or the 
expected steady-state profit denoted by g may be obtained after allowing t 
tend to infinity from (2.5.15). 

Therefore 


(^ 1 +-“1)0^2 +/^ 2 ) 


(2.5.16) 


In order to examine economic behaviour of the system in a finite 
interval of time t (t is not large enough) let us specify various parameters as 
below : 


^3 = 0.050, = 0-025, = 1.500, = 1.000, 

roo = 175.0, r„ = 50.0, r 22 = 75.0, r 33 = -215.0, 

ry = 0 for i j. 

Substituting these values in (2.5.15) we obtain ¥o(t) for various values 
of t as given in table 2.5.1 : 


70 



t 

Voft) 

t 

Voft) 

0.0 

0.0 

8.0 

1352.33147 

0.5 

86.61606 

9.0 

1520.73082 

1.0 

172.07458 

10.0 

1689.12988 

2.0 

341.49547 

11.0 

1857.52884 

3.0 

510.19627 

12.0 

2025.92776 

4.0 

678.68960 

13.0 

2194.32666 

5.0 

847.11942 

14.0 

2362.72556 

6.0 

1015.52878 

15.0 

2531.12446 

7.0 

1183.93130 




Table 2.5.1 


From above table we see that for t > 12 the increase in value of Vo(t) 
for every unit increase in value of t is 168.39890 which is nothing but the 
steady-state expected profit of the system. 

Profit Optimization Further if one is posed with the problem of 
selecting failure rate and repair rate jx from a given set of these parameters 
along with cost parameters, we should proceed as follows. 

(a) Specify cost matrix (rjj), values of ^ and jU corresponding to each 
alternative. 

(b) Develop cost equations for the system model under consideration 
and obtain expression for g. 

(c) Choose that alternative taking values of /S and fi as decision 
variables for which g is maximum. 

In what follows, we present the results in table 2.5.2 for the system 
under consideration. 
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Table 2.5.2 


From the above table we observe that for the sixth pair of units 
steady-state expected profit is maximum. 

Particular Case 

We now discuss a particular case of the model when both the units 
are identical. Let the failure and repair rates be ^ and fi respectively. Define 
the following three states 

State Condition of the system 

50 Both the units are operating. 

51 One unit operates, the other is in failed state and is under repair. 

52 Both the units are in failed state and are under repair. 
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The corresponding cost matrix for such a system can be given by 

^00 ^01 ^02 

^10 ^11 ^12 

^20 ^21 ^22 


Obviously r^j = 130 = 0. Since these transitions are not possible, results 
for this particular case can easily be obtained for the model discussed above 
after making the following observation. 


Original model 

Particular case 

^1- h 

P 

H 

P- 

’’ov '’02 

r01 

’’lO' '’20 

'’10 

^31' ''32 

^21 

'’i3' '’23 

^12 

r-[i, ^22 

•’ll 

'■33 

^22 



Therefore 


g=777TTl^^(’'oo+^roj)+^/f(rji+//rio+^rj2)+^V22'''^^2i)] (2.5.17) 

(P+/“) 

To see the analytic behaviour of g in (2.5.17), let us consider the 
following cost structure : 


foo ~ 


% = fi(«) = 

^22 ~ ^2 ^2 


A > 0 

aj > 0, Cj < 0 

0 , (^ < 0 


and 







So 


‘01 


^lo = r 


12 


‘21 


0. 


8 “ (2-5-18) 

{p-^n) 

To determine that vaue of fi which maximizes expected profit, g in 
(2.5.18) for a given failure rate and known values of cost parameters, we 
differentiate g with respect to fi and put equal to zero, i.e., dg/cl^ = 0, 
which gives the value of denoted by jii which maximizes expected profit, 
g denoted by g* 

2(A+2^c,)-(2a,+^C2) 


since is finite and cannot be negative, being repair rate, to decide fi* gives 


d^e 

maximum we have to further examine the sign of 


(i) 


14. ={4 


If / 2 * > 0 and 


ft =f4 


, < 0, is a maximum point for 


g i.e., the sufficient conditin for /i* > 0 to be absolute maximum 
point for g is 

( A + ^Cj ) “ (2<z j +j8c2) ^ 


• d^B 

(ii) However, if w >0 and — 

dpr 


fi =fi 


. > 0, /f* is the minimum 


point for g i.e., when (A+^Cj)-( 2 o:j+^C 2 )+a 2>0 


and absolute maximum is achieved at f4*=Q for which g =0^. 
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CHAPTER HI 

OPTIMAL FAILURE ANALYSIS IN 
SERIES SYSTEMS 



3.1 INTRODUCTION 


The maintenance problems related to complex systems with many 
component modules are quite difficult. Analysis of failures plays an important 
role in the evaluation studies of maintained systems. Keeping in view the 
complexity and economics of modern systems, it is of paramount interest to 
spend on different phases (parts) of a system in an optimal manner. In 
industrial/defence systems, whenever a failure occurs, it is normal to first 
decide what type of failure it is. More elaborately, specialization and diversity 
in repair/or nature of failure demand certain type of activities to be performed 
during the time lag, when a failure occurs and the instant when it really 
enters repair. These activities may be termed as ‘failure analysis’. A close 
look at their behaviour leads to the following main activities to be covered 
under failure analysis. 

(i) Detection of failure This may be the case with systems which 
malfunction or operate with reduced efficiency on the occurence 
of a failure and the problem is to detect a failure/fault. 

The concept of fault detection in maintained systems is important 
and has been considered by workers in different contexts. In 
fact, faults may be of 3 types ; 

(a) Identifiable faults are those faults which are identified as and 
when they occur. In other words, sometimes the subsystems responsible for 
malfunctioning are observable and failure diagnosis in such cases is straight 
forward. 

(b) Unidentifiable faults are those faults which caimot be identified 
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without any outside intervention. These are hidden faults and may cause 
serious damages in the long run. Such faults can only be identified/detected 
by the help of detectors. However, detectors may be used in the case of 
identifiable faults as well. In other words, there are systems which shut down 
completely as soon as failure occurs in one of its constituent component 
modules and such systems do not show any discernible symptoms for detecting 
the faulty component module. A serially connected heater serves a good 
example. 

(c) Recoverable faults are those faults which are automatically 
recovered, once they occur. This feature of maintained systems is taken care 
of at design stage itself. 

(ii) Identification and location of failure Once the failure is 
detected the next problem which arises is that of failure 
identification and its location i.e., to identify the phase or module 
which has failed. 

(iii) Classification of failure Having arrived at the defective 

component, one may be interested to classify the failure in one 
of the several prespecified classes. These classes may depend 
upon the severity of failure, time/cost to be spent on repair. 

(iv) Administrative time on failure rectification There may be 
some administrative problems related to failure/repair. The time 
and money spent on these activities are also a part of failure 
analysis. 

(v) Transportation of the faulty part for repair > The resources 
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spent on transporting the faulty module (part) to repair shop 
also constitute a part of failure analysis. 

In order to have a good understanding of the genuine requirements of 
the above activities, let us consider the example of a dockyard in which ships 
arrive from time to time for refit. The normal procedure adopted is when a 
ship reports the dockyard for repair, a reviewing comittee meets and decides 
about the type of repairs, time/cost estimates needed on the ship at various 
work centres in the dockyard. Later on, the ship goes to different work centres 
and repair is carried out. Here, the time spent on administrative affairs and 
other activities pertaining to failure analysis is quite substantial and need be 
included. 

The above activities and several other not listed above, need attention 
of a system analyst responsible for system analysis/evaluation studies of the 
system under consideration. 

A number of authors have concenterated their attention to reliability 
parameters after including some of the above activities. Kumar and Aggarwal 
(1978) have analyzed a 2-unit standby system with negligible time spent on 
fault analysis along with two types of repairs. They obtained MTSF, steady-state 
availability, absorption probabilities, expected number of visits to a state, 
second moment of time spent in an up state etc. Kapoor and Kapur (1975b) 
obtained pointwise availability for two models of a complex system with two 
classes of components and 1 ^. They assumed a probability distribution for 
the time spent on post repair. Lai et al. (1980) analyzed a 2-unit standby 
system with ‘analysis time’. This ‘analysis time’ is constituted either of failure 



analysis time or of post repair time. The two models have been analyzed in 
detail and several parameters of interest are obtained. 

There are not may articles on failure-analysis from the view point of 
cost analysis. Nevertheless, there are articles dealing with optimization problems 
after taking due consideration of underlying costs. 

Gross (1970) considered a 2-unit series systems with failure detection 
mechanism. He associated a probability of correct detection of faulty component 
which depends upon the cost spent on the detection mechanism. It is assumed 
that in case of wrong detection of failure (i.e., assigning failure to a unit 
which has not failed in reality) there is some extra cost spent on the repair 
of a good unit. He minimized the expected misclassification cost and obtained 
a condition under which a detection mechanism would be economically feasible. 
Later, Kumar (1975) extended the results to an n-unit series system. It is 
interesting to note that probabilities of detection of failure have been taken 
as functions of cost. The properties of these cost functions have already been 
cited in Chapter 1. 

Takami et al. (1980) have discussed the problem of allocation of fault 
detectors to an n-unit series system. The detectors are subject to failure viz., 
they produce unwanted alarms and they do not produce alarms when they 
are needed. System equations have been written and steady-state availability 
has been obtained. Further, assuming cost of detectors and loss due to system 
failure, optimum number of detectors has been determined to minimize 
expected cost by applying 0-1 programming. Kumar and Kapoor (1979) 
analyzed a 2-unit series model with detectors which are subject to 2 types of 
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failures. 


All the papers referred above deal with fault detection. As regards 
fault classification, there has not been much work. Kumar and Kapoor (1981b) 
have developed a classification rule based on life-time distribution for a 2-unit 
series system which minimizes expected cost of misclassification. Butterworth 
(1972) has considered fault testing modules for complex systems assuming a 
priori pobabilities of modules. Park (1982) has developed a mathematical 
model for diagonising a faulty component module in an n-unit series system. 
He has suggested an optimal testing policy based on cost incurred on testing 
and causative failure probabilities of the underlying modules. 

In the present chapter we concentrate on the following important 
problems pertaining to failure analysis in maintained systems. 

1. We have developed a decision rule based on life time distribution 
of the system that minimizes the expected misclassification cost 
in an n-unit series system. 

2. A 2-unit series system with 2 detectors has been examined with 
a view to allocate cost to detectors which maximizes the expected 
steady-state profit of the system. 

We discuss the above in detail in the subsequent sections : 

3.2 OPTIMAL FAILURE CLASSIFICATION IN AN n-UNIT 

SERIES SYSTEM 

We discuss below the problem of classification of failures in a series 
system. 
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Failure Classification System Model and Problem Definition 

1- There is a series system consisting of n units. Let the units be 
denoted by 1, 2, n. The units are independent and have 
general distributions for their failure times. Obviously, the system 
fails upon failure of any one of the n-units. 

2. At the time of system failure, the problem is to classify this 
failure in one of the n classes (mutually exclusive), each class 
corresponds to a unit in the system. Thus, when a failure is 
classified into i‘^ class it implies that the system failure has 
occured due to failure of the 1“’ unit. While dealing with failure 
classification, one is prone to the error of classifying a failure. 
More elaborately, a failure ought to be classified in class i, is 
classified in class j. This leads to loss, may be termed as 
‘misclassification cost’. 

3. Further, a priori probabilities of failure of units are assumed to 
be known. 

Notation 

i, j : index, which denote units label = 1, 2, ..., n. 

Cjj : cost of misclassification of a failure, assigning the failure to 

unit j whereas unit i has failed. 

fi(t) : life-time pdf of unit i. 

R : sample space 

Rj : subset of R which corresponds to failure of unit i. 


80 



R = U R. 

i=l 

Py(R) : probability of misclassification of a failure due to unit j 
whereas unit i has failed. 

= S f,(t)dt 

R 

j 

: a priori probability of failure of unit i. 


Optimal Classification : Loss Minimization and Classillcation Rule 

Development 

The expected loss due to misclassification is given by 

E(L) = i [ i C,jP,(R)] (3.2.1) 


The problem is to choose RjS’ so that E(L) is minimized. Since we 
know a priori failure probabilities for each unit and corresponding life-time 
pdfs, we may define the conditional probability of failure due to unit i as 



k=l 


(3.2.2) 


On classifying the failure due to unit j, the expected loss is 


i^j i Vk(t) 


C. 


(3.2.3) 


This expected loss is to be minimized by the proper choice of j. 
Equivalently, the problem is to minimize 
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(3.2.4) 



The choice of j may be made by evaluating (3.2.4) for all values of j’s 
and choose the one which gives the minimum. If there are more than one 
j’s giving the same minimum, any one of them may be chosen. In this manner 
RjS are formed. Thus the classification procedure is to classify failure to unit 
j if it falls in R. 

J J 

In what follows is the following theorem : 


Theorem If Jt- is the a priori probability of a failure due to unit i 
in an n-unit series system with probability density function fj(t) and if is 
the cost of misclassification assigning the failure to unit j whereas unit i has 
failed then the regions of classification, Rj, Rj, ...., R^, that minimizes the 
expected loss are defined assigning the failure to k belonging to R^ if 


2 C„ < 2 

1 = 1 1=1 

i5^k 


^ifi(t) C-- , = 


(3.2.5) 


(If (3.2.5) holds for all j (j^^k) except for h indices and the inequality 
is replaced by equality for those indices, then this point can be assigned to 
any one of h + 1 units). 

Discriminatory Power of the Classification Rule Now we suggest a 
method to study the discriminatory power of the classification rule developed 
above. 

Let ajj denote the number of failures known which have occured due 
to failure of unit i but which were classified as have occured due to failure 
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of unit j. Then the matrix for the classification problem is defined to be 
m X m matrix S = {a-^. this matrix may be called the RightAVrong Identification 
Matrix (RWIM) and is given by 



■«n 


^13 

• • * ^Im ' 



^21 

^22 

^^23 

®2ni 


s = 

• 

• 

• 

• 

(3.2.6) 


An.1 

«m2 

«m3 

• • • 



Diagonal elements of S denote the number of correct classifications 
and off diagonal elements denote the number of misclassifications. 

Then, we define the normalised RWIM, Sq as 

So = 

^ ij m 

I «ij 

Evidently, the elements of normalised RWIM give proportions of 
classifications : right and wrong. 

For the purpose of testing the discriminatory power of the classification 
procedure, we make use of a test based on distribution. Accordingly, we 
define 

Q = tars)! + (3.2.7) 

e e 

where n and n denote the number of right and wrong classifications made 
by the discriminatory procedure; e and e denote the expected number of right 
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and wrong classifications that would be made if the classifications were made 

m 

at random. Then, if N = 2 N. denotes the total number of classifications 

1=1 

made, and the probability of a successful random classification is = , then 

m 


m 

3 

II 

n = 

N-n 


N 


.T N 

e = — , 

e = 

N- — 

m 


m 


On substituting (3.2.8) in (3.2.7), we get 

O = (N-nm)^ 

^ N(m-l) 


(3.2.8) 


(3.2.9) 


a form convenient for numerical evaluation. Now, if we define : 

Hq ; classification is made at random 

Hj : classification is made making use of (3.2.5). 

Under Hg, (3.2.9) follows Chi-square distribution with (m^-1) degrees 
of freedom. 


Hq will be rejected or accepted accordingly as calculated value of 
Chi-square, is larger or smaller than the tabulated value. 

In order to explain the usage of the results developed above, we now 
give a numerical example. 


Illustration 

Let us consider a 3-unit series system with units 1, 2 and 3 and with 
exponential life-time distributions for the units, i.e., 
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fi(t) = i = 1, 2. 3 

Then (3.2.5) reduces to 


I c„ < i q 

1 = 1 issl 

i^k 

We assign failure to unit 1 if the following set of inequalities are satisfied 

^2^2(^)^21 ^3^3(0^31 < ^1^1 (0^12'^ ^3^3(0^32 

^2^2(^)^21 "*'^3^3(^)^31 ^/i(0Ci3+^2^2W^23 

The failure is assigned to unit 2 if the following set of inequalities are 

satisfied. 

^1^i(^)^ 12’^^3^3(^)^32 ^lfl(l^)Cj3+^2^2(^)^23 

^l^l(^)^12'^^3^3(^)^32 ^2^2(^)^21 '^^3^3(^)^31 

The failure is assigned to unit 3 if the following set of inequalities are 

satisfied. 



i ^ ^ 2 ^ 2 ( 0^21 "^'^ 3 ^ 3 ( 0^31 

^2^2(^)^23“^^1^i(^)^13 

; ^lfl(l)Ci2'*"^3^3(^)^32 

Now for given set of values 


= 0.185 

C 21 = 20 . 0 , 

C 33 = 30.0 

^2 = 0.040 

C 22 = 15.0, 

C 32 = 25.0 

^3 = 0.080 

Ci 3 = 15.0, 

C 23 = 20.0 

We find JTj by 
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(i = l,2,3) 


(3.2.10) 



i = l 


On substituting values of and in (3.2.10), we obtain 

TTj = 0.60656, 7t^ = 0.13115, = 0.26230 

Also making appropriate substitutions in inequalities (3.2.5), we see 
that when t > 40.2360 component 2 is classified as failed, when 
9.3670 < t < 40.2360 component 3 is classified as failed and when 
o < X < 9.3670 component 1 is classified as failed. 


Now to test the discriminatory power we say, let the 100 classifications 


are made which are given in RWIM as below 



1 

2 

3 

Totals 

1 

10 

5 

7 

22 

2 

6 ! 

25 

9 

40 

3 

16 

10 

12 

38 


Then, we observe N = 100, n = 47, m = 3 

(lOO^MZ)! = 8 405 
^ 100(3-1) 


Also Chi-square at 1% level is 1.646 is certainly significant. Tlius the 
classification procedure certainly does better than chance. 


33 OPTIMAL ALLOCATION OF COST TO DETECTORS IN A 2-UNIT 
SERIES SYSTEM 
System Model 

(i) There is a 2-unit series system. Let the units be identified as 1 
and 2. 
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(ii) Each unit is equipped with a separate detector to detect failure. 
As soon as the system failure occurs i.e., one of 2-units fails, 
detectors will detect the failures with some prespecified 
probabilities. 

(iii) There is no interconnection between the detectors. In fact, 
probability of correct detection is directly related to the cost 
spent on detectors i.e., more the cost spent on the detector, the 
better probability of correct detection (smaller probability of 
wrong detection) will be ensured. Some of the properties of 
these probabilities are discussed in Kumar (1975). 

(iv) When a unit is wrongly detected as failed, it enters for repair. 
It is assumed here that the repair is not carried out completely 
but after a short (random) duration it is identified that there 
was wrong detection. As such the failed unit goes to repair 
without any ambiguity now. 

(v) When the system is down, the operative unit cannot fail. 

(vi) It is assumed that detectors do not fail at the time of need. 

(vii) Failure time distributions of units are assumed to be exponential 
whereas repair time distributions are general. 

(viii) The system earns (loses) at a fixed rate in each state which can 
be different for each state. There is a fixed transition reward 
(cost) whenever the system changes its state. 

Define following system states at any instant. 
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System states and transitions 

Sq : Both the units 1 and 2 are in operative state. 

Sj : Unit 2 fails but unit 1 is detected as failed and goes for repair. 

$2 : Unit 1 fails but unit 2 is detected as failed and goes for repair. 

53 : Unit 2 fails and is detected as failed and goes for repair. 

54 : Unit 1 fails and is detected as failed and goes for repair. 

The system is up in Sq and is down in Sj (i = l, 2 , 3 , 4 ). Transitions 
between various states are given in figure 3 . 3 . 1 . 



Figure 3 . 3.1 Transitions between various states 


Notation 

ySj : constant failure rate of unit 1 

^^2 : constant failure rate of unit 2 
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nij 

m2 


Pnih) ■ 

P2l(^2) • 

Pllih) ■ 

Pllih) ■ 

i i : 
Pij = 

p=[pij) 


mean repair time of unit 1 
mean repair time of unit 2 

mean time for which the repair is carried out on unit 1 in 
Sj wrongly whereas unit 2 was actually in failed state or mean 
time taken for detecting unit 2 as failed. 

mean time for which the repair is carried out on unit 2 in 
Sj wrongly whereas unit 1 was actually in failed state or mean 
time taken for detecting unit 1 as failed. 

probability of detecting unit 2 as failed whereas unit 1 has 
failed when Xj units of money are spent on detector attached 
to unit 1. 

probability of detecting unit 1 as failed whereas unit 2 has 
failed when X 2 units of money are spent on detector attached 
to unit 2. 

probability of detecting unit 1 as failed when unit 1 fails and 
Xj units of money are spent on detector attached to unit 1. 

probability of detecting unit 2 as failed when unit 2 fails and 
X 2 units of money are spent on detector attached to unit 2. 

subscripts that imply system states i, j = 0, 1, 2, 3, 4. 

One step transition probability for making a transition from 
S, to 

transition probability matrix 
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I : identity matrix of order 5 

D : I - P 

dj : i"’ subdeterminant of D, deleting i"’ row and i”’ column 

TTj : steady-state probability in Sj = d/2 dj 

n : row vector, » ^ 4 ) 

Tj : mean unconditional waiting time in Sj 

Tjj : earning rate in S,- 

Tjj : transition reward (cost) for making a transition from Sj to Sj 

g : total expected earning per unit time in steady-state or 

steady-state expected profit. 

Expected Profit of the System 

Observing state transitions from figure 3.3.1, various state transition 


probabilities, pjj can 

be written 

as 


Poi 

00 

= P2l(X2)J' 

0 

y32e~^2*e“^>‘ 

dt=P2l(X2) 



P02 

00 

= Pl2(Xl){ 


dt=Pi2(Xi) 

A 

(3.3.1) 

P03 

00 

= P22(X2)/ 

0 

)32e~^2‘e~^i‘ 

dt=p22(X2) 

^1+^2 


Po 4 

00 

= Pll(xl)/ 


dt=Pii(Xi) 

A 

A+A 



Obviously, pjj = P24 = P30 = P40 " ^ 
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Further, it is easy to see that 


To= 7 dt= 

0 

Tj = nijj, T2 = 


T3=m2, 


Also 


T^ = Tn^ 


^21 ^22(^) ~ ^ 

Following Barlow and Proschan (1965), 
we know 


71 = 2rP 

2 r(I-P) = 0 

so, we obtain dj as follows 

do=l, dj=poi, d2=po2, d3=pQi+po3, d4 = Po2 + Po4 

Also 


(3.3.2) 


(3.3.3) 


(3.3.4) 

(3.3.5) 


(3.3.6) 


by 


di 



d; 


(3.3.7) 


j 

Following Howard (1971), steady-state expected profit, g can be given 


* “ '2 Ji.T, 


(3.3.8) 


where 
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Substituting (3.3.7) in (3.3.8), we get 

I diT,q, 

® “ 2 <1,T, 


On making appropriate substitutions in (3.3.9) and taking rij = 0 for 
also using (3.3.3) and after simplification g is obtained as 

l+^2[P2l(X2H2 + n^2]+A[Pl2(Xi)m„ + mj] 1 • • 7 

If no cost is spent on detectors, it is reasonable to think that on the 
failure of the system, any unit (1 or 2) can go for repair with equal probability 
i.e.. 


Pi2(0) =P2i(0) =Pii(0) = P22(0) =1 


Then 


^2 

roo+^2r33ni2 +^1^44^1 +Y rumi2+Y r22mii 

IXItO 1^71 

1 


(3.3.11) 


(3.3.12) 


Optimal Cost Allocation to Detectors 

Now, if one is posed with the problem of optimal cost allocation to 
detectors attached to units 1 and 2 with a view to maximize expected profit, 
g of the system under given cost constraint, then the problem reduces to 
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Maximize g in (3.3.10) 
subject to 

Xj + Xj = X Xj, X 2 , X > 0 (3.3.13) 

where Xj, Xj are the costs to be spent on detectors attached with units 1 and 
2 respectively, x is the total amount available. 

We now discuss the above problem in the light of two cost structures 
used for failure detection mechanism. 


Failure Detection Cost Structure 1 
Let 


P12 (Xi) 


1 

2(1 +Xj) 


and 


P2l(X2) = 


1 

2(1+X2) 


(3.3.14) 


It may be noted that (3.3.14) satisfies properties given in Kumar (1975). 

On substituting (3.3.14) in (3.3.13), the problem is to maximize g subject 
to Xj+X 2 =x where g is given by 


roo+^2r33m2+^ir44”^i+2(i+x2) 2(1+X2) 

1 +^1^1 + 2( 1 +X2) ^ 2(1 +Xj) 


Observing the equality, x^ + Xj = x, (3.3.15) reduces to 


(3.3.15) 
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(3.3.16) 


Axj+I3xj+E 

Bxj+Fxj+G 


where 


and 


A = -2Ai, D = 2AjX + Dj-Ej, E = (l+x) (2Aj+Ei) + Dj 

B = -2Bi, F = 2BjX + Fi-Gj, G = (l +x) (2Bj+Gi) + Fj 

■^1 ”■ ^00 "^^ 2 ^ 33 ^ 2 2 ^11 ^ 12 ’ ~^ 1^22 

Bj = 1 +^2^12 2 ^12’ Gj =^jnijj 

To maximize g in (3.3.16) we differentiate g with respect to Xj and put 


(Jp 

^=0, which gives 


-L±V L^-KM 
K 


(3.3.17) 


where 


K=AF - DB, L = AG - EB, M = GD - EF. 

As Xj is cost we consider only positive value of Xj, and denote it by 


To examine if x/ is maximum point, it is ethical to verify 


dx? 


Xi=Xi 


<0 


Once Xj* is determined, (the optimal cost spent on the second 
* * 

detector) is obtained as X 2 = x - Xj . 
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Illustration 


To illustrate the results numerically, let us specify the various parameters 


for the problem as under : 




o 

11 

/52 = 0.5, 

3 

11 

o 

bo 

mj = 0.4, 

mjj = 0.1. 

mj2 = 0.2, 

^00 ~ 20.0, 

r„ = -22.0, 

^22 ~ "23.0, 

r33 = -25.0. 

r44 = -27.0. 

x = 50.0. 




When the above values 

are substituted 

in (3.3.17), we 

get 


Xj* = 15.242, Xj* = 34.758 
and 

g = 7.8083 

On the other hand, if no cost is spent on detectors, (3.3.12) gives 
g = 6.5845 

Thus inclusion of detectors ensures an increase of 18.58% in expected 

profit. 


Failure Detection Cost Structure 2 
Let 

Pi2(xi)=| 

and 

P2i(x2)=| 


(3.3.19) 
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where aj and aj are constants depending upon the requirement of correct 
detection assurance. 


On substituting (3.3.19) in (3.3.13), the problem is to maximize g subject 
to Xj + X 2 = X where g is given by 


^2 - 

roo+^2h3^2+^i^u^i+-2^ “"^rijm^+ye “'**r22mij 
1 +^e"‘‘2^mi2+-:^e'®f'im,i 


(3.3.20) 


Observing the equality Xj + X2 = x, (3.3.20) reduces to 


g 


A+Be~Vi+He^^^ 

E+Fe-®i".+Ke"^i 


where 


^00 ^2^33^2 ^1^44^1’ 


B = 


^ 1 ^ 22^11 T' ^ 1^11 -ir ^ 2 ^ 11 ^ 12 ® ^ T, ^ 2 ^ 12 ® ^ 


, F 


, H 


E = 1 +^ 2 ^! + 
, K 


(3.3.21) 


The value of a2X has been bounded for the traceability of results i.e., 

I a2X i < 4. 

To maximize g in (3.3.21) we differentiate g vwth respect to Xj and put 
dfi 

^ = 0, which gives 

Sie'^A+S2e^+S3e<"2'"i^i=0 (3.3.22) 

where 

Si=ai(AF-EB), S2 = a2(EH-AK), S3 = (aj + a2) (FH-KB). 
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Let 


t = and ai = 32 in (3.3.22) 

Then 

S2t^ + S3t + Si = 0 (3.3.23) 

Solving (3.3.23) for t, we get 

-S3±V s^s^s" 

2S2 

Since Xj is cost, we consider only that value of t which gives Xj > 0 

from 



where log is taken to the base e. 

Let us denote this value of Xj by Xj*. 

To decide that Xj* actually maximizes g in (3.3.21), it is ethical to verify 

. < 0 
Xi=X, 

Once Xj’ is determined, Xj* (the optimal cost to be spent on second 
detector) is obtained as X 2 = x - Xj’. 

Illustration 

To illustrate the results numerically, let us specify the various parameters 
as under : 
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= 0 . 02 , 


m, =0.8, 



g = 184.3361. 

On the other hand if no cost is spent on detectors, we find from 
(3.3.12) 

g = 182.2168 

Thus, we observe that employing detectors increases the expected profit 
of the system. 
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CHAPTER IV 


EVALUATION OF SERIES 
SYSTEMS THROUGH DYNAMIC 
PROGRAMMING TECHNIQUE 



4.1 INTRODUCTION 


Designing reliability into a system is often required when the system 
is to perform an important mission for industry/military or space flight and 
where malfunction (failure) of a constituent part causes heavy losses in terms 
of money/time. 

In order that a multistage system having many subsystems or stages in 
series have a failure free operation, the subsystems’ reliability should be very 
high. While dealing with evaluation studies of series systems, one often finds 
following ways to increases system reliability. 

(a) To increase the reliability of the constituent subsystems 
(components) i.e., to develop highly reliable components. This 
in turn increases cost of the system. Further, as the number of 
series components increases, the overall system reliability 
decreases. In fact, after a certain stage, even a marginal increase 
in a component’s reliability may result in tremendous cost. 

(b) To allocate certain components in redundancy in an optimal 
manner at various stages of a series system so that the overall 
system reliability is increased. 

(c) To suggest some maintenance policies for an existing system with 
a view to increase system reliability. 

A little concentration on the above discussion reveals that (a) and (b) 
are useful at the design stage, whereas (c) plays an important role in the 
evaluation/maintenance problems of the system. More elaborately, the methods 
developed to meet the requirements of (a) and (b) lose their importance, 
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once a system is designed and put into operation. Tlius, at the latter stage 
viz., operation stage, one must look for methods that may help in better 
maintenance of the system. Besides, problems at the maintenance stage are 
still more important because, sometimes state of the art does not allow a 
further increase in components’ reliability inspite of one’s willingness to bear 
extraordinary amount of money. 

There has been a large number of papers in the above directions. 
[Barlow and Proschan (1965), Bellman and Dreyfus (1958, 1962), Kettle (1962), 
Proschan and Bray (1965), Fyffe et al. (1968), Hadley (1964), Messigner and 
Shooman (1970), Jensen (1970), Woodhouse (1972) and Misra (1971a). 

Misra (1971a) gave a dynamic programming formulation to determine 
optimal number of redundant components at each stage in a series system. 
Lagrangian Multiplier Technique was used to solve the problem. Initially, the 
problem was formulated under one constraint. However, redundancy allocation 
under two constraints has also been dealt with. 

Bellman and Dreyfus (1962) have discussed the reliability R of a system 
containing several subsystems, each of which has a redundancy. They maximized 
the reliability subject to a cost constraint. Kettle (1962) has discussed an 
approach of least cost allocation with a view to minimize sum of subsystems 
unreliabilities instead of maximizing system reliability. This was extended by 
Proschan and Bray (1965) and Barlow and Proschan (1965) to the case of 
two constraints. 

In formulating and solving redundancy allocation problems through 
dynamic programming with a view to optimize system reliability, number of 
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constraints causes difficulty because of increased dimensionality. This difficulty 
has been successfully overcome by employing Lagrangian Multiplier Technique 
by Bellman and Dreyfus (1958), Fyffe et al. (1968), Hadley (1964), Messigner 
and Shooman (1970) . Kettle (1962) and Woodhouse (1972) have used the 
concept of dominating sequences for the purpose. 

However, (c) has also attracted attention of some of the workers. Kumar 
and Lai (1980) studied a 2-unit cold standby redundant system and determined 
an optimal maintenance policy which maximizes expected profit of the system 
assuming that there is no discounting. Later on Kumar and Kapoor (1981a) 
extended the procedure to cover the case of discounted cost structure. 

It may be worthwhile to point out that in most of the papers cited 
above, the methods of analysis are based on dynamic programming formulation. 
In fact, dynamic programming is a very powerful technique and its spectrum 
is quite wide in that it enables solution procedures for a fairly large reliability 
situations. Howard’s (1971) policy iteration method is also a version of dynamic 
programming. 

In the present chapter we discuss the following problem for series 
systems, viz.. Optimum cost allocation under a constraint on the cost in the 
subsequent section. 



4.2 OPTIMUM COST ALLOCATION IN A SERIES SYSTEM UNDER 
A CONSTRAINT ON THE COST 
System Model 

There is an n-component series system. The problem is to allocate cost 
to each component with a view to maximize system reliability subject to total 
cost. It is assumed that there exist a number of components with vaiying costs 
and reliability which are capable of performing the same system task. The 
reliability of each component is assumed to be monotonically increasing 
function of its cost. 

Notation 

R„(C„): reliability of component n 
Cp : cost of component n 

R : system reliability 

N : number of components in the system. 


(1961) has 

(4.2.1) 


System Reliability 

Then, the system reliability R is given as 
R(c.. C, C„) = n R„(C.) 

n=l 

Now, we obtain a suitable expression for Rn(C„). Beriphol 
expressed the cost of a component as a function of its reliability 
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where and Kjn are constants. 

Equation (4.2.1) cannot be solved for R (C ). 

n' ft' 

Therefore, we devise a function given below in (4.2.2) which leads to 
approximately the same set of values for reliability and cost as is given by 
(4.2.1) for 0.9 < R„(CJ < 1. 

Rn(C„) = a„ + b„logC„ (4.2.2) 

where 0 < a„ < 1, b„ > 0 and log is to the base e. 

Now the actual problem becomes 

Maximize R(Cj, C 2 ,...., C*.;) 

^ 2 ’ > 

subject to 

i C„=C , C„ > 0 

n=l 

N 

i.e., Maximize [ 0 (a„+b„ log C„) j 

» 

subject to 

i C„ = C, C„>0 (4.2.3) 

n=l 

Dynamic Programming Formulation 

With a view to maximize system reliability in N stages, define stage j 
(j = l, 2, ..., N) to consist of j components counted from the last i.e., N to 
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(N-j + 1). 

Define 

F(j, Cf^_j+i, C) = Reliability of stage j when cost is allocated 

to component N-j + 1, is optimally allocated to the remaining j-1 

components contained in stage j, j = 1, 2,..., N. 

F (j, C) = Maximum of F (j, C) with respect to 

Cn-j+i. 0 < <C, j = 1, 2, N. (4.2.4) 

Evidently, F* (N, C) = Max [R(Ci, C 2 ,..., Cj„)] 

with respect to Cj, C 2 ,...., Cj,. 

subject to 

+ C2 + + Cf^T — c 

Bellman’s Principle of optimality results in the following functional 
equations : 

F* 0',C) = Max F* (j-1, C-C^^^,)] 

with respect to , 0 < Cj^_j+i<C, j = 2, N 

F’ (1, C) = Max [R^,(C^)] 

with respect to Cjs^ , 0 < < C (4.2.5) 

The optimum values of Cj, C 2 , ...., obtained by solving above 

functional equations. 

If standard dynamic programming procedures are followed to solve 
(4.2,5) the solution is quite voluminous and cumbersome. Therefore, it is 
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desirable to achieve approximate solutions. We describe below a simple 
approach to obtain approximate but explicit expression for optimal costs 
through an illustration. 

Illustration 

Let 

N = 3 , 

From ( 4 . 2 . 5 ) 

F ( 1 , C) = Max [R3(C3)] with respect to C3 , 0 < Cj < C 

= Max [aj + b3 log C3] with respect to C3 , 0 < C3 < C 
= 83 + b3 log C ( 4 . 2 . 6 a) 

Evidently, 

F( 2 ,C 3 , C) = R3(C2) F* ( 1 , C-Q 

= i^+h , log q) [a3 + b3 log (C-q)] ( 4 . 2 . 6 b) 

Let 

K = C^C ( 0 <K< 1 ) 

gi = aj + bj log C, i = 1 , 2, 3 
Then 

F( 2 , CK, C) = g2g3 + b3g2 log ( 1 -K) + b2 g3 log K + b2b3log K. log ( 1 -K) 
Neglecting bjbj log K log ( 1 -K) as 62 and bg are small quantities. 

F( 2 , CK, C) = g2g3 + b3g2 log ( 1 -K) + b2g3 log K = f(K) 

The mflYiTnnm value of K denoted by K* is obtained by differentiating 
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^p(K) with respect to K and equating it to zero i.e., = 0 which gives 

Qlx 


K 


^283 


^382'*'^2S3 


To decide that K provides maximum, we further examine the sign of 


dV(K) 

dK^ 


K=K* 


i.e., 


d > (K ) 

dK^ 

Thus, 

dV(K) 

dK^ 


|.^2& ^3§2 1 

(1-Kr 


K=K 


. < 0 


So K* maximizes ^(K) and hence 
F*(2, O = <p (K*) = F(2, CK‘, C) 
On similar lines for third stage 


F(3, Cj, C) = Ra(Q)F* (2, C-C,) 

= [aj+bj log CJ [a 2 + b 2 log (C-Cj) + bj log K ] . 
[a 3 +b 3 log (C-Cj) + b 3 log (1-K*)] 


Let 


L=Ci/C (0 < L < 1) 

Then 
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where 


F( 3 , CL, C) = [aj + bi Jog C + bj Jog L] [a2 + b2 log C + bj log ( 1 -L) 
+ bj log K ] [ aj + bj log C +b3 log ( 1 -L) 

+ b3 log ( 1 -K*)] 

= [hj+bi log L] [h2 + b2 log ( 1 -L)] [ h3 + b3 log ( 1 -L)] 


= gi 

^2 = g2 + ^>2 log ^ 

1^3 = g3 + ^3 log (1-K‘). 

So, on neglecting the terms containing bjb2, b2b3, bib3 and bjb2b3 as 
bj, b2, b3 are small quantities. 

F( 3 , CL, C) = hih2h3 + hjh2b3 log ( 1 -L) + hib2h3 log ( 1 -L) 

+ bj h2h3 log L 
= m (say). 

The maximum value of L denoted by L* is obtained by differentiating 
^(L) with respect to L and equating it to zero i.e., 

Ma = 0 

dL 

which gives 

_ l^il^2l^3 

bih2h3+b2hjh3 + b3hjh2 

To decide L’ provides maximum, we further examine the sign of 
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dL' L=L 


* i.e., 


d^g(L) _ 


b^hjha-Fbghih; 

V (1-Lf 


Thus . < 0 

dL^ L=L* 


So, L maximizes ^(L). 

Optimum values of Cj, C 2 , C 3 denoted by Cj*, C^r , are obtained 


as follows 


Cj = CL, C 2 '= CK* (1-L’), C 3 ' = C(1-K') (1-L‘). 


Illustration 

We now illustrate the results numerically. Let us take C = 15,00,0(X). 
Further, assume the following data taken from Beriphol (1961) 


Ri 

Cl 


^2 

R, i 

3 j 

C3 

.999 

.995 

1820 

100 

.999 

.995 

1030 

100 

.999 

.995 

6310 

200 


Using the above set of values of and C„ (n = 1 , 2, 3), the values 
of a„ and b„ (n = l, 2, 3) are obtained from (4.2.2) as 

aj = .97910, 32 = .97530, 83 = .98086 

bj = .00138, b 2 = .00171, b 3 = .00116 

Consequently, 

K’ = .5937, L* = .3246 
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and 


Cl* = 4,86,900; q* = 6,03,503.67; C 3 * = 4,09,596.33 

Evidently in this case the maximum attainable reliability of the system 
due to optimum allocation of cost to each component for a given total cost 
of 15 lakhs is .99111. 
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